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Abstract 

This  thesis  examines  a  remote  sensing  technique  for  measuring  the  atmospheric  struc¬ 
ture  constant  (C^)  and  the  transverse  atmospheric  wind  velocity  (V)  as  a  function  of  alti¬ 
tude  by  performing  temporal  and  spatial  correlations  of  band  pass  filtered  wave  front  slope 
sensor  measurements.  Two  point  sources  are  used  to  illuminate  two  pairs  of  co-located  wave 
front  slope  sensors.  The  sources  and  sensors  are  arranged  to  give  rise  to  crossed  optical 
paths.  For  each  pair  of  wave  front  slope  sensors,  the  output  of  the  larger  sensor  is  sub¬ 
tracted  from  the  output  of  the  smaller  sensor.  This  band  pass  filters  the  smaller  sensor’s 
output.  The  filtered  outputs  are  then  correlated.  It  is  shown  that  the  spatial  correlation  is 
related  to  the  structure  constant  by  an  integral  containing  and  a  path  weighting  func¬ 
tion.  The  path  weighting  function  is  sharply  peaked  at  the  intersection  altitude  and  decays 
rapidly  to  zero.  It  is  also  shown  that  the  temporal  correlation  is  related  to  V  through  an 
integral  containing  V,  C^,  and  a  wind  path  weighting  function  which  is  sharply  peaked  at 
the  intersection  altitude  and  decays  rapidly  to  zero.  The  vertical  resolution  of  each  path 
weighting  function  is  calculated.  Also,  signal-to-noise  ratios  are  calculated  using  generally 
accepted  models  for  the  structure  constant  and  atmospheric  wind  profiles. 
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REMOTE  SENSING  OF  TURBULENCE  AND  TRANSVERSE  ATMOSPHERIC 


WIND  PROFILES  USING  OPTICAL  REFERENCE  SOURCES 

/.  Introduction 

1.1  Atmospheric  Turbulence  and  Winds 

Lightwaves  traveling  from  space  to  Earth  are  distorted  when  they  pass  through  the 
Earth’s  atmosphere.  This  distortion  gives  rise  to  the  well  known  twinkling  (properly  termed 
scintillation)  of  stars,  and  is  the  limiting  factor  in  astronomers’  efforts  to  obtain  high 
resolution  images  of  celestial  objects.  In  theory,  the  resolution  achievable  by  an  imaging 
system,  called  the  diffraction  limited  resolution,  is  proportional  to  the  size  of  the  system’s 
aperture  (5).  Yet  atmospheric  distortion  reduces  the  resolving  power  of  large  aperture 
telescopes  to  that  of  a  telescope  with  a  10  to  20  cm  aperture  (14). 

A  growing  field  of  study  called  adaptive  optics  attempts  to  restore  the  telescope’s 
theoretical  resolution  by  compensating  for  the  atmospheric  distortions.  Generally,  adaptive 
optics  systems  incorporate  deformable  mirrors  to  remove  wave  front  distortions  (8).  A 
control  system  senses  the  incoming  wave  front  distortion  and  sends  commands  to  the 
mirror.  The  mirror  assumes  a  shape  that  is  the  conjugate  of  the  wave  front  distortions. 
Thus,  distortions  are  removed  when  the  wavefront  reflects  off  of  the  mirror. 

The  efficient  design  and  operation  of  adaptive  optics  systems  require  characterization 
of  the  distribution  of  turbulence  along  the  wave  front’s  propagation  path.  This  distribution 
can  be  determined  from  a  parameter  called  the  atmospheric  structure  constant.  The 
structure  constant  is  discussed  further  in  the  next  section. 

Along  similar  lines,  meteorologists  and  aviators  are  interested  in  knowing  the  atmo¬ 
spheric  wind  velocity  as  a  function  of  altitude.  Such  information  helps  increase  weather 
forecasting  accuracy  and  aviation  safety  (13).  Additionally,  knowledge  of  the  vertical  wind 
profile  is  used  in  the  determination  of  control  loop  bandwidth  for  adaptive  optics  systems. 
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In  both  cases,  measuring  the  winds  at  specific  locations  within  the  atmosphere  is  impor¬ 
tant  and  can  be  accomplished  by  determining  the  wind  profile  along  an  optical  path  which 
passes  through  the  points  of  interest. 

1.2  Turbulence  and  Wind  Characterization 

Differential  heating  of  the  Earth  causes  atmospheric  turbulence.  This  heating  gener 
ates  pockets  of  air,  called  eddies,  with  varying  temperatures  and  pressures,  giving  rise  to 
varying  indexes  of  refraction.  The  index  of  refraction  is  denoted  by  n,  and  is  a  func.'on  of 
both  time  and  position  within  the  atmosphere.  Thus,  the  fluctuation  oi  n  is  a  temporal 
and  spatial  random  process  (25).  Generally,  the  temporal  fluctuations  are  assumed  to  be 
due  to  atmospheric  winds,  an  assumption  known  as  Taylor’s  hypothesis  (19).  Taylor’s 
hypothesis  essentially  states  that  the  pockets  of  air  do  not  decay  as  they  are  blown  across 
an  optical  system’s  line  of  sight. 

The  strength  of  the  turbulence  is  characterized  by  the  spatial  power  spectral  density 
(psd)  of  the  fluctuations  of  the  index  of  refraction.  The  psd  is  denoted  by  $„(/v'),  where 
K  is  the  spatial  wavenumber  vector.  The  magnitude  of  ^n{K)  is  proportional  to  the 
atmospheric  structure  constant,  C^.  In  many  optical  propagation  geometries,  the  structure 
constant  is  a  function  of  position  along  the  optical  path.  This  is  particularly  true  for  vertical 
propagation  through  the  atmosphere.  Due  to  its  relationship  with  the  psd,  is  used  as 
a  measure  of  the  strength  of  the  turbulence  along  the  optical  path. 

As  alluded  to  above,  the  temporal  fluctuations  of  n  are  assumed  to  be  due  to  atmo¬ 
spheric  winds.  Wind  has  direction  as  well  as  magnitude,  and  therefore  is  denoted  by  the 
vector  V .  If  the  time  it  takes  for  an  eddy  of  known  size  to  blow  across  the  optical  path 
can  be  measured,  then  the  wind  velocity  can  be  determined.  Note,  however,  that  only  the 

component  of  the  wind  perpentlicular  to  the  optical  path  can  be  measured  in  th,.,  _ niier. 

As  will  *^e  discussed  later,  the  temporal  fluctuations  of  n  can  be  used  to  derive  information 
about  the  tn  nsverse  wind  velocity  as  a  function  of  position  along  the  optical  path. 
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1.3  Measurement  of  and  V 

No  instrument  can  measure  directly.  Likewise,  although  anenometers  can  measure 
V,  such  devices  cannot  be  placed  at  all  the  points  necessary  to  obtain  a  good  profile  of 
V  along  an  optical  path.  The  usual  remote  sensing  approach  is  to  measure  some  effect  of 
the  atmosphere  on  a  reference  beam  of  light  and  relate  the  measured  effect  to  or  V. 
Numerous  techniques  have  been  used  to  measure  and  V  (1,  3,  10,  21).  Such  techniques 
have  become  more  practical  with  recent  advances  in  the  generation  of  artificial  guide  stars 
for  use  as  reference  sources  (4).  Most  measurement  techniques  have  exploited  the  spatial 
or  temporal  cross-correlation  properties  of  the  intensity  of  the  optical  field.  The  research 
presented  here  takes  a  different  approach,  one  based  upon  the  cross-correlation  properties 
of  the  optical  field’s  phase. 

1.4  Approach 

The  basic  measurement  geometry  to  be  examined  is  shown  in  Figure  1.1.  Two  wave 
front  sensors  separated  by  AS  are  present  in  the  aperture  plane  of  an  optical  system.  The 
wave  front  sensors  are  illuminated  by  two  point  sources,  pi  and  p2,  which  are  separated 
by  Ap.  The  geometry  is  arranged  to  give  rise  to  crossed  optical  paths.  The  wave  front 
distortions  on  the  two  reference  beams  caused  by  the  turbulence  at  the  intersection  point 
will  be  highly  correlated,  while  distortions  caused  by  turbulence  elsewhere  along  the  paths 
will  not  be  correlated.  This  fact  is  used  to  determine  C^.  The  intersection  altitude  can  be 
varied  by  changing  the  source  or  sensor  separations,  and  in  this  way  a  vertical  profile  of 
the  turbulence  strength  can  be  obtained. 

Additionally,  time-lagged  correlations  of  the  wave  front  distortions  will  be  well  corre¬ 
lated  when  the  lag  equals  the  time  necessary  for  an  eddy  to  traverse  the  distance  separating 
the  two  beams.  Thus,  the  time-lagged  correlations  of  the  wave  front  distortions  on  the  ref¬ 
erence  beams  can  be  used  to  determine  V. 

Rejack  examined  the  measurement  geometry  depicted  in  Figure  1.1  for  use  in  mea¬ 
suring  C"’  (17).  He  demonstrated  that  is  related  to  the  correlation  of  the  two  wave 
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Altitude  z 


Figure  1.1.  Basic  Measurement  Geometry  (17) 


front  sensor  measurements  by 


{C,}  =  {s,it)s,{t)) 

=  j  dz' C^iz')wciz')  (1.1) 

where  (•)  represents  the  ensemble  average  operator,  Sn{t)  is  the  output  of  wave  front  sensor 
n  at  time  t,  z'  is  the  position  along  the  z  axis,  and  Wc(z')  is  called  the  path  weighting 
function.  Rejack  then  improved  the  system’s  resolution  performance  by  linearly  combining 
measurements  from  two  sets  of  wave  front  sensors. 

This  thesis  investigates  the  crossed  beam  measurement  geometry  of  Figure  1.1,  but 
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examines  a  new  method  for  combining  the  measurements  from  the  two  sets  of  wave  front 
sensors.  There  are  three  basic  objectives  of  this  research.  The  first  objective  is  to  de¬ 
rive  the  new  method’s  path  weighting  function  wd^)  and  evaluate  its  characteristics. 
The  second  objective  is  to  extend  the  new  technique  to  the  remote  sensing  of  transverse 
atmospheric  winds.  It  is  shown  that  the  derivative  of  the  time-lagged  correlation  of  the 
measurements  can  be  related  to  V(z)  by  (21) 

^(C,(r))|^=0  =  ■^{Siit)S2{t  +  T))\r:=0 

=  j  dz'Cl{z')V{z')-wviz')  (1.2) 

where  ^{C,)  is  the  partial  derivative  with  respect  to  r  of  the  time  lagged  correlation  of 
the  sensor  measurements,  is  the  band  pass  filtered  output  of  wave  front  slope  sensor 
n,  T  is  the  amount  of  time  lag,  V{z')  is  the  wind  velocity  at  the  altitude  z',  and  Wv(z') 
is  called  the  wind  path  weighting  function.  The  wind  path  weighting  function  is  derived 
and  its  characteristics  examined.  Finally,  the  third  objective  is  to  develop  signal-to-noise 
ratio  expressions  which  will  help  determine  how  many  measurements  are  needed  to  make 
an  accurate  estimate  of  the  quantity  of  interest.  Following  the  fulfillment  of  these  three  ob¬ 
jectives,  an  example  is  presented  to  demonstrate  the  and  wind  path  weighting  function 
shapes  obtained  when  a  practical  measurement  geometry  is  used.  Finally,  some  practical 
implementation  issues  and  their  impact  on  system  performance  are  discussed. 

1.5  Overview 

Chapter  II  contains  a  review  of  atmospheric  turbulence  and  wind  sensing  techniques. 
Chapter  III  contains  the  methodology  used  in  conducting  this  research,  including  limiting 
assumptions.  Chapter  IV  consists  of  a  signal-to-noise  ratio  analysis  of  the  measurement 
technique  proposed  in  this  thesis.  Chapter  V  presents  an  example  measurement  geome¬ 
try  and  discusses  some  implementation  issues  which  affect  system  performance.  Finally, 
chapter  VI  provides  conclusions  and  some  suggestions  for  further  research. 
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II.  Background 


2.1  Theory  and  Definitions 

2.1.1  Turbulence  Atmospheric  turbulence  is  the  result  of  differential  heating  of  the 
Earth.  This  differential  heating  causes  the  formation  of  various  sized  pockets  of  air,  called 
eddies.  The  atmosphere  consists  of  a  random  distribution  of  these  eddies,  each  with  a 
different  temperature  and  pressure.  The  index  of  refraction,  n,  is  defined  as  the  ratio 
of  the  speed  of  light  within  a  medium  to  the  speed  of  light  within  a  vacuum  (7).  The 
speed  of  light  within  the  atmosphere,  and  therefore  the  index  of  refraction,  is  a  function  of 
temperature,  humidity,  and  pressure.  Thus,  as  different  portions  of  an  optical  wave  front 
pass  through  different  eddies,  the  wave  front  sections  travel  at  different  speeds.  This  speed 
differential  distorts  the  propagating  wave  front. 

2.1.2  Power  Spectral  Density  of  n  The  fluctuation  of  the  index  of  refraction  within 

the  atmosphere  is  characterized  by  a  wide  sense  stationary  random  process  and  there¬ 
fore  has  an  associated  power  spectral  density  (psd),  denoted  as  $„(A^),  where  K  is  the 
wavenumber  vector  and  is  a  measure  of  the  number  of  eddies  of  size  L  =  (6).  For 

isotropic  turbulence,  that  is,  turbulence  for  which  the  autocorrelation  of  the  fluctuation  of 
n  has  spherical  symmetry,  K  is  replaced  by  if  =  \K\. 

In  a  classic  work,  Kolmogorov  divided  ^n(A')  into  three  separate  regions  based  upon 
eddy  size  (9).  The  first  region  contains  eddies  of  size  greater  than  Lq,  where  Lo  is  called  the 
outer  scale  size.  The  shape  of  ^„iK)  in  this  region  is  not  predictable  by  his  theory.  The 
second  region,  known  as  the  inertial  subrange,  contains  eddies  between  the  sizes  Lq  and 
io,  where  £o  is  called  the  inner  scale  size.  The  form  of  $„(A')  within  the  inertial  subrange 
can  be  predicted  by  the  theory  of  turbulent  flow,  and  is  given  as 

$„(A')  =  0.033C^K^  (2.1) 

where  C\  is  the  structure  constant.  The  third  and  final  region  contains  eddies  smaller 
than  lo-  The  psd  in  this  region  decays  rapidly  towards  zero.  When  modeling  atmospheric 
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turbulence  effects  on  optical  systems,  the  inertial  subrange  is  generally  the  only  region  of 
concern. 

2.1.3  Index  of  Refraction  Structure  Function  The  structure  function  of  n  is  com¬ 
monly  used  to  characterize  the  second  order  moment  of  n.  The  structure  function  is  defined 
as  (6) 

r^)  =  {[n{rl)  -  n(r;)]^)  (2.2) 

where  ri  and  are  position  vectors  and  (•)  is  the  expected  value  operator.  If  the  turbulence 
is  homogeneous,  meaning  it  is  spatially  wide  sense  stationary,  as  well  as  isotropic,  and  if  it 
has  a  psd  given  by  the  Kolmogorov  expression  of  equation  2.1,  then  the  structure  function 
becomes  (6) 

D„{r)  =  Clri  (2.3) 

where  r  =  |rj  -  rl|. 


2.2  Turbulence  and  Wind  Profile  Models 

Now  that  the  atmospheric  structure  constant,  C^,  has  been  related  to  turbulence 
strength  through  use  of  the  structure  function,  some  information  about  is  needed.  The 
atmospheric  structure  constant  cannot  be  measured  directly.  However,  experimental  data 
has  been  employed  to  derive  and  validate  some  models  for  the  vertical  profile  of  C^.  One 
of  the  models  most  often  used  is  called  the  Hufnagel- Valley  turbulence  model  and  has  the 
form  (16) 


Cl{z')  =  5.94  X  10“®^  f — )  2'“’  exp  f 

’  \2l)  *^Vl000/ 


(2.4) 

where  z'  is  altitude  in  meters,  w  is  the  average  wind  speed  in  m/sec,  and  A  is  a  free 
parameter  adjusted  to  yield  the  desired  profile  shape.  The  units  of  are  m"^.  Equation 
2.4  is  plotted  in  Figure  2.1  for  w  =  27  and  A  =  10“’'*  if  z'  >  11000,  0  if  ^'  <  11000. 


Unlike  C^,  the  wind  velocity  can  be  measured  directly,  and  models  of  the  wind 
velocity  profile  have  been  created  based  upon  empirical  data.  The  wind  velocity  profile 
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Figure  2.1.  Hufnagel- Valley  Profile  (16) 


most  often  used  in  turbulence  problems  is  the  Bufton  model  given  by 


V(z‘)  =  Vg  30  exp 


-  9400 \ 
4800  ) 


(2.5) 


where  z'  is  altitude  in  meters,  and  Vg  is  ground  wind  speed  in  m/sec.  Equation  2.5  is 
plotted  in  Figure  2.2  for  Vg  =  5  m/sec. 


0  5000  10000  15000  20000 

Altitude  (Meters) 

Figure  2.2.  Bufton  Wind  Profile  (16) 
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2.3  Evaluation  of  and  V 

The  technique  most  often  used  to  remotely  sense  and  V  is  to  extract  these  values 
from  the  correlations  of  optical  scintillations.  The  scintillation  data  is  obtained  by  measur¬ 
ing  the  irradiance  intensity  produced  by  a  source  located  within  or  above  the  atmosphere. 
Fried,  in  1969,  was  one  of  the  first  to  examine  the  relationship  between  the  spatial  covari¬ 
ance  of  intensity  and  Cl  (3).  By  including  temporal  effects.  Fried  was  also  able  to  relate 
the  temporal  covariance  of  the  optical  scintillations  to  V .  Because  his  work  is  the  basis  for 
many  similar  measurement  techniques,  Fried’s  approach  is  summarized  in  the  next  section. 

2.3.1  Correlations  of  Irradiance  Measurements  Fried  proposed  measuring  the  in¬ 
tensities  incident  on  two  telescopes  using  the  geometry  shown  in  Figure  2.3.  The  telescopes 
have  highly  restricted  fields  of  view  and  receive  a  monochromatic  signal  from  a  stellar 
source.  The  intensities  measured  at  time  t  are  denoted  ii{t)  and  *2(0-  The  logarithmic 
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Distorted  Wavefront 

\ 

\ 
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Telescope  1 


Telescope  2 


Figure  2.3.  Fried’s  Measurement  Geometry  (3) 


amplitude,  Imit),  is  calculated  using 


Ut)  =  ^\n 


.(0 


(UJ 


(2.6) 
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where  (i^)  is  the  average  Vcilue  of  The  spatial-temporal  logarithmic  amplitude 

covariance  is  then  given  by 


t)  =  ([/i(0  -  (^>][^2(i  +  7-)  -  (/2)])  (2.7) 


where  li{t)  and  /2(<)  are  the  log-amplitudes  as  given  in  equation  2.6,  and  p  is  the  separation 
of  the  two  telescopes. 

Fried  suppresses  the  time  dependence  by  setting  r  =  0  and  applies  Tatarski’s  result 
that  (18) 

C/(p)  =  0.652A:^  y  dz  C^{z)  J  d<7  (7~  Jo(o’p)  ^1  -  cos  (2.8) 

where  k  =  z  is  the  distance  from  the  telescope  to  a  point  along  the  propagation  path,  L 
is  the  distance  from  the  receiver  to  the  source  (in  this  case  infinity),  and  Jq  is  the  zero-order 
Bessel  function. 

Through  a  series  of  integrations,  Fried  reduces  equation  2.8  to 

CM  =  ki  jf”  iz  z>CHz)T  (^)  (2.9) 

where  (^)  groups  all  of  the  terms  that  depend  upon  (^)-  Fried  tabulates  T 
for  a  few  values.  Thus,  Fried  has  been  successful  in  relating  to  the  measured  irradiance 
correlations  through  a  linear  integral  equation. 

Fried  next  includes  the  time  dependence  of  the  intensity  fluctuations  to  relate  the 
temporal  covariance  of  the  log-amplitudes  to  atmospheric  wind  velocities.  Thus,  equation 
2.8  becomes 


Ci(p,  t)  =  0.652k^  J  dzCl(z)J  d<7  Jo((7\p  -  V  (z)t\)  1 1 


a^z 


cos 


(2.10) 


where  V{z}  is  the  wind  velocity  at  altitude  z. 
Accordingly,  equation  2.9  becomes 


CMr)  =  kif  izz>Ci(z)r(^^}EzJ^ 


(2.11) 


2-5 


As  can  be  seen,  knowledge  of  both  Ci{p,T)  and  Cl(z)  is  necessary  to  solve  for  K(^). 
Even  then,  determining  V{z)  requires  solving  a  nonlinear  integral  equation,  a  difficult 
if  not  impossible  task.  However,  Fried  states  that  numerical  techniques  will  yield  an 
approximate  solution,  although  the  procedure  requires  an  initial  estimate  of  V(z).  An 
interesting  and  important  point  is  that  no  data  beyond  that  already  obtained  to  determine 
Cn{z)  is  necessary  to  determine  V{z}.  However,  additional  processing  is  required. 

2.3.2  Crossed  Beam  Technique  In  1974,  Wang  et  al.  (21)  used  a  crossed  beam 
technique  to  profile  both  C^(z]  and  V{z).  In  1991,  Beland  and  Krause-Polstorff  examined 
the  same  crossed  beam  technique,  but  used  range  gated  lidar  pulses  to  generate  optical 
sources  within  the  atmosphere,  allowing  vertical  profiling  of  C\  (1).  Wang’s  analysis  is 
based  upon  the  phase  screen  approach  of  Lee  and  Harp  (11).  The  measurement  geometry 
studied  by  Wang  et  al.  is  shown  in  Figure  2.4. 


The  measurement  geometry  was  formed  by  two  laser  sources,  Tq  and  Tj,  separated  by 
p(,  and  two  receivers,  Ra  and  flj,  separated  by  xq,  j/o-  The  transmitters  and  receivers  were 
arranged  in  such  a  way  that  the  optical  paths  from  Ta  to  Rg  and  from  Tj  to  intersected. 
The  2  axis  was  perpendicular  to  the  planes  containing  the  transmitters  and  receivers. 
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Given  this  measurement  geometry,  Wang  et  al.  found  that  the  amplitude  covariance 
function  is  given  by 

--  J  dz'C^{z')W,i,iz')  (2.12) 

where  Cab  is  th®  covariance  of  the  intensities  measured  by  and  Rb,  and  Wab  is  called 
the  C\  path  weighting  function  and  is  given  by 


Woi(2')  =  0.1327r“A:'''  /  dA'A'  sin 

Jo  L  J  L  I  ^  \ 

(2.13) 

where  L  is  the  distance  from  the  transmitter  plane  to  the  receiver  plane,  and  z'  is  the 
distance  from  the  transmitter  plane  to  some  point  along  the  z  axis. 


U<^s{L  -z')^ 

Jr. 

K 

2kL  \ 

Jo 

(xo  -  yo)z'  ,  ^ 

- ? - +  Pi  -  T 


Wang  and  his  colleagues  found  a  similar  expression  relating  the  amplitude  covariance 
function  to  the  transverse  wind  velocity.  They  accomplished  this  by  inserting  the  time 
dependence  into  the  intensity  covariance,  and  then  defining  a  new  function,  Mabi  to  bo 


^ab  — 


dCab{puT) 


dr 


T=0 


=  /  dz'CUz')Viz')-W:„iz') 
Jo 


(2.14) 


where  Wa^{z')  is  called  the  wind  path  weighting  function  and  is  expressed  as 


/•OO 

/  dA'A'^  sin^ 

rA'V(^-^')l 

/o 

2kL 

Ji(Kr)f 


(2.15) 


The  function  is  the  first  order  Bessel  function,  and  f  is  a  unit  vector  along  f  where  f  is 
expressed  as 

f  =  (zo  -  jTo)^ (2.16) 

An  important  aspect  of  this  technique  is  that  it  does  not  measure  the  wind  component 
parallel  to  the  z  axis.  This  is  because  an  eddy  moving  parallel  to  the  z  axis  would  intercept 
both  beams  at  the  same  time,  thus  having  no  effect  on  the  slope  of  the  amplitude  covariance 
at  zero  time  lag. 
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2.3.3  Correlations  of  Phase  Measurements  Rather  than  measuring  intensity  cor¬ 
relations,  Welsh  (22)  has  proposed  using  the  spatial  correlation  properties  of  the  optical 
field’s  phase  to  profile  C^(z).  Welsh  examines  a  crossed  beam  measurement  geometry  us¬ 
ing  artificial  guide  stars  as  the  atmospheric  sources  and  wave  front  slope  sensors  as  the 
receivers.  This  geometry  is  shown  in  Figure  2.5.  Welsh  shows  that  the  correlation  of  the 


z 


Figure  2.5.  Measurement  Geometry  Examined  by  Welsh  (22) 


slope  sensor  outputs  is  given  by 

(C*)=  rdz'Cliz')wciz')  (2.17) 

Jo 
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where  the  weighting  function  Wc{z')  depends  upon  the  source  altitude,  z,,  source  separa¬ 
tion,  Pi  —  p2,  receiver  separation.  Ax,  PSD  of  the  refractive  index,  and  the  sensor  aperture 
weighting  functions.  He  also  includes  a  term  to  account  for  any  linear  post  processing  that 
may  be  performed  on  the  wave  front  sensor  data. 

Welsh,  after  assuming  a  Kolmogorov  spectrum  for  evaluated  the  path  weighting 
function  numerically  for  two  cases.  First,  the  weighting  function  is  evaluated  for  the  case 
correlating  the  direct  phase  measurements  of  the  wave  front  sensors.  In  the  second  case,  the 
wave  front  sensor  phase  measurements  are  band  pass  filtered  before  being  correlated.  Welsh 
found  that  the  band  pass  filtering  provided  a  more  sharply  peaked  weighting  function, 
which  translates  into  better  altitude  resolution.  Such  a  result  is  not  surprising.  The  low 
frequency  phase  effects  remain  well  correlated  further  from  the  intersection  point  because 
these  effects  are  caused  by  large  eddies.  Band  pass  filtering  removes  these  low  frequencies 
and  thus  helps  the  phase  correlation  function  to  decay  more  quickly,  improving  altitude 
resolution. 

Rejack  (17)  further  expanded  upon  Welsh’s  technique.  Whereas  Welsh  assumed  the 
wave  front  sensors  actually  measured  the  wave  front  phase.  Rejack  considered  a  more 
realistic  wave  front  sensor  in  which  the  overall  wave  front  slope  is  measured.  Rejack 
began  by  deriving  an  expression  for  the  path  weighting  function  wc(z)  that  depended 
only  upon  the  measurement  geometry.  In  accomplishing  this,  he  assumed  a  Kolmogorov 
spectrum.  Rejack  plotted  Wciz)  and  showed  that,  like  Welsh’s  no  post  processing  case, 
Wc(z)  possessed  a  .slow  decay  rate  yielding  poor  altitude  resolution. 

Rather  than  post  process  the  wave  front  sensors’  data.  Rejack  investigated  a  tech¬ 
nique  for  combining  the  weighting  functions  resulting  from  the  correlations  of  slope  mea¬ 
surements  from  two  different  size  apertures.  Specifically,  he  subtracted  the  phase  corre¬ 
lation  function  for  a  pair  of  apertures  of  size  L2  from  that  for  a  pair  of  apertures  of  size 
Li,  where  Li  <  L^-  Rejack  plotted  the  resultant  modified  weighting  function  for  three 
different  ratios  of  Li  to  L^.  He  found  that  the  new  weighting  function  did  decay  more 
rapidly  than  that  due  to  a  single  pair  of  apertures. 
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2-4  Conclusion 

Remote  sensing  of  atmospheric  parameters  using  optical  beams  is  a  widely  studied 
‘ield.  The  crossed  beam  technique  has  been  used  to  obtain  correlations  of  optical  fields 
and  thus  derive  information  about  the  distributions  of  and  V.  While  many  of  these 
methods  have  relied  upon  correlations  of  intensity  scintillations,  recent  works  have  also 
examined  the  use  of  phase  correlations  in  the  crossed  beam  technique  for  sensing  C^.  It  is 
proposed  that  the  time  lagged  phase  correlation  properties  of  the  crossed  beam  technique 
may  be  used  to  obtain  a  vertical  profile  of  V(z').  It  is  also  proposed  to  combine  the 
phase  measurements  of  co-located  wave  front  sensor  pairs  in  a  manner  which  improves  the 
resolution  of  the  measurement  technique  beyond  that  previously  achieved.  The  following 
chapter  contains  the  development  of  these  proposals. 
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III.  Methodolofpj 


3. 1  Introduction 

As  discussed  in  previous  chapters,  most  technirpies  for  sensing  C'“  and  V  have  relied 
upon  measurements  of  optical  scintillations.  Some  work  has  also  explored  the  use  of  mea¬ 
surements  of  wave  front  phase  to  sense  C^.  This  chapter  presents  a  method  for  combining 
wave  front  phase  measurements  from  co-located  wave  front  phase  slope  sensors  to  improve 
the  altitude  resolution  achievable  by  the  crossed  beam  measurement  geometry.  Addition¬ 
ally,  this  chapter  shows  that  wave  front  phase  measurements  may  also  be  used  to  obtain 
vertical  profiles  of  V. 

Figure  3.1  illustrates  the  measurement  geometry  considered  in  this  chapter.  Two 
point  sources,  denoted  as  pi  and  p2  and  separated  by  Ap,  are  present  at  altitude  z,.  These 
point  sources  illuminate  four  wave  front  sensors  which  are  present  in  the  aperture  plane  of 
an  optical  system.  The  sensors  are  denoted  Si,  Sj/,  s^,  and  So'.  Slope  sensors  and 
are  centered  on  the  same  point  in  the  aperture  plane.  The  vector  distance  from  the  center 
of  one  sensor  pair  to  the  center  of  the  other  sense  pair  is  Ax.  The  field  of  view  of  each 
wave  front  sensor  is  restricted  so  that  sensors  Sj  and  Si<  receive  light  from  pi  only,  and 
sensors  S2  and  S2'  receive  light  from  p2  only.  As  can  be  seen  from  Figure  3.1,  this  geometry 
gives  rise  to  crossed  optical  beams.  The  altitude  of  the  crossing  point  is  controlled  by  the 
source  and  sensor  separations. 

In  the  measurement  technique  examined  here,  the  output  of  wave  front  sensor  is 
subtracted  from  the  output  of  s„.  This  is  done  for  both  n  =  \  and  n  =  2.  The  subtraction 
essentially  band  pass  filters  the  slope  measurement  of  which  Welsh  (22)  has  shown 
to  be  a  necessary  step  in  obtaining  a  sharply  peakf  !  weighting  function.  The  filtered 
measurements  are  then  correlated.  This  correlation  is  designated  (C,)  and  is  defined  by 

{C,)  =  ((-Sl  —  5l')(^2  ~  -52'))  (^-l) 

The  goal  of  the  first  derivation  presented  in  this  chapter  is  to  relate  (C,)  to  C^.  For  this  first 
derivation,  the  beams  physically  intersect  and  all  measurements  are  taken  simultaneously. 


3-1 


Altitude  z 


Figure  3.1.  Me«isurement  Geometry  (17) 

The  first  derivation  results  in  an  equation  of  the  form 

{C,)  =  j  Az'Cl{z')wc{z')  (3.2) 

where  z'  is  a  position  along  the  z  axis  and  vjc(z')  is  called  the  path  weighting  function. 
Ideally,  Wc(z')  would  have  the  shape  of  a  Dirac  delta  function  so  that  the  sifting  property 
could  be  employed  to  relate  the  measured  correlation  directly  to  C^. 

In  the  second  derivation  presented,  the  measurements  are  not  made  simultaneously. 
Instead,  the  measurements  on  sensors  S2  and  33'  are  taken  r  seconds  after  the  measurements 
on  sensors  Sj  and  Sj'.  Additionally,  the  optical  beams  are  oriented  so  that  there  is  always 
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a  separation  between  the  two,  even  at  the  cross  over  altitude.  Thus,  there  is  a  small 
separation  of  the  beams  at  the  crossing  altitude  and  the  beams  do  not  physically  intersect. 

Given  these  two  conditions,  the  second  derivation  shows  that  the  derivative  of  the 
time  lagged  correlation  of  the  wave  front  slope  measurements  is  related  to  V  by  (21) 

^(Cs(T)}jr=0  =  -  Sl'(t)][s2(f  -t)-  S2<(i  -  t)]) 

=  Jdz'Cl{z')Viz')-nry{z')  (3.3) 

where  ^{C,{r))  is  the  partial  derivative  with  respect  to  r  of  the  time  lagged  correlations 
of  the  sensor  measurements,  r  is  the  amount  of  time  lag,  V(z')  is  the  wind  velocity  at  the 
altitude  z\  and  wv{z')  is  called  the  wind  path  weighting  function. 

3.2  Justification  of  Approach 

As  mentioned  above,  the  outputs  of  co-located  slope  sensors  are  subtracted  before 
any  correlations  take  place.  This  subtraction  has  the  effect  of  band  pass  filtering  the 
slope  measurement  of  the  smaller  sensor.  Welsh  (22)  has  demonstrated  that  this  approach 
provides  a  more  sharply  peaked  path  weighting  function  than  that  found  when  no  filtering 
is  performed. 

To  understand  how  this  band  pass  filtering  is  accomplished,  it  is  necessary  to  know 
something  about  the  frequency  response  of  wave  front  slope  sensors.  An  idealized  wave 
front  slope  sensor  frequency  response  is  shown  in  Figure  3.2.  The  frequency  response  has 
an  upper  cutoff  frequency  inversely  proportional  to  the  size  of  the  slope  sensor  (22).  Thus, 
the  smaller  the  wave  front  slope  sensor,  the  higher  the  cutoff  frequency.  From  this,  it  is 
easy  to  see  how  subtracting  the  measurements  of  a  larger  slope  sensor  from  those  of  a 
smaller  co-located  slope  sensor  will  result  in  a  band  pass  filtering  of  the  signal  from  the 
smaller  sensor.  The  low  frequencies  common  to  the  measurements  of  both  sensors  will 
be  eliminated,  while  those  frequencies  existing  between  the  cutoff  frequency  of  the  larger 
sensor  and  the  cutoff  frequency  of  the  smaller  sensor  will  be  unaffected.  Figure  3.2  presents 
a  graphical  depiction  of  this  technique. 
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Figure  3.2.  Band  Pass  Filter  Creation 


In  Rejack’s  analysis  (17),  he  correlated  measurements  from  the  sensors  Si  and  Sq, 
and  then  correlated  measurements  from  the  sensors  Si/  and  52-.  He  then  investigated  the 
difference  of  the  two  correlations 

(C,}  =  (S1S2}  —  (si'S2')  (3.4) 

This  approach,  however,  is  not  equivalent  to  the  band  pass  filtering  method  suggested 
by  Welsh.  In  fact,  Welsh  (22)  showed  that  such  an  approach  is  inferior  to  the  approach 
adopted  in  this  work  where 

(C,>  =  ((si  -  Si-)(S2  -  S2'))  (3.5) 

As  has  been  discussed  previously,  obtaining  a  vertical  profile  of  the  wind  velocity 
requires  time  lagged  correlations  of  the  phase  slope  measurements.  The  time  lagged  cor¬ 
relation  is  similar  to  equation  3.5,  but  with  the  addition  of  a  time  lag  for  measurements  S2 
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a:id  S2' 

(^.(r))  =  {[si(0  -  Si-(<)][«2(«  -  r)  -  S2>(t  -  r)])  (3.6) 

where  r  is  the  time  delay  between  measurements. 

S.3  Derivation  of  Path  Weighting  Function,  wc{z') 

3.3.1  Overview  The  first  objective  of  this  research  is  to  find  an  expression  for  the 
path  weighting  function  which  results  when  the  slope  sensor  outputs  are  combined  as 

shown  in  equation  3.5.  Once  the  form  of  this  new  weighting  function  is  determined,  the 
performance  of  the  measurement  technique  can  be  examined.  Accordingly,  the  goal  of  this 
derivation  is  to  relate  the  correlation  shown  in  equation  3.5  to  through 

{€,)  =  I  dz'C^Jz')wciz')  (3.7) 

where  z'  is  a  position  along  the  2  axis  and  Wciz'}  is  called  the  path  weighting  function. 
In  accomplishing  this  goal,  the  first  step  is  to  write  an  expression  for  the  correlation  defined 
in  equation  3.5.  This  requires  use  of  a  model  for  the  output  of  the  wave  front  slope  sensors. 
Then  Kolmogorov  statistics  are  assumed  for  the  atmospheric  turbulence.  This  allows  the 
path  weighting  function  to  be  reduced  to  an  expression  which  depends  only  upon  the 
measurement  geometry.  Thus,  it  is  then  possible  to  examine  the  characteristics  of  wc(z')- 

3.3.2  Slope  Signal  Correlation  First  consider  the  single  realization  correlation  de¬ 
fined  as 

C's  =  (Si  —  .Si/)(52  —  52')  (3.8) 

where  s„  is  the  output  of  slope  sensor  n.  Slope  sensors  Sj  and  Si-  are  co-located  and  are 
separated  from  52  and  52'  by  Ax.  The  subscripted  numbers  correspond  to  the  illuminating 
light  source,  either  pi  or  p2-  Normally,  wave  front  slope  sensors  consist  of  multiple  sub¬ 
apertures.  However,  it  will  be  assumed  throughout  this  development  that  the  wave  front 
slope  sensors  are  each  composed  of  just  one  aperture. 

Let  (l>„{x)  designate  the  phase  of  the  incident  optical  field  from  source  p„.  The  output 
of  wave  front  slope  sensor  n  is  a  noisy  measurement  of  the  average  slope  of  d>„(x)  across 
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the  sensor’s  aperture  and  is  denoted  s„.  The  expression  fc>r  s„  is 

Sn  =  j  [V<^„(f)  •  d]  +  a„  (3.9) 

where  V0„(x)  is  the  spatial  gradient  of  the  phase,  Wn{x)  is  the  sensor’s  aperture  weighting 
function,  d  is  a  unit  vector  in  the  wave  front  sensor’s  direction  of  sensitivity,  and  a„  is 
the  slope  measurement  error.  The  measurement  error  is  caused  by  photon  noise  in  the 
detector.  The  photon  noise  is  sometimes  referred  to  as  shot  noise. 

Wallner  integrated  equation  3.9  by  parts  to  yield  (20) 

Sn  =  -  j  d^x  W’(x)(f>„(x)  +  a„  (3.10) 

where  IT^(x)  is  the  gradient  of  Wn{x)  in  the  wave  front  sensor’s  direction  of  sensitivity. 

Referring  to  Figure  3.1,  it  will  be  assumed  that  the  inner  sensors,  described  by  Wi 
and  VV2,  are  identical  except  for  being  placed  at  different  positions.  The  same  assumption 
is  applied  to  Wii  and  IV2'.  Using  these  assumptions,  the  notation  can  be  simplified  by 
letting  VFi(x)  =  W2{x)  —  Wi^{x),  and  Wi<(x)  =  W2'{x)  =  Wi^{x)  from  here  on.  The 
subscripts  Li  and  L2  refer  to  the  dimensions  of  the  apertures. 

With  the  aid  of  the  slope  sensor  output  model  given  by  equation  3.10,  equation  3.5 
may  be  written  as 

Cs  =  [[-J  d^x  W’,(f)d>i(x)  +  ai  -  [-J  dV  W’,(x')(/>j(x')  +  Qi.]| 

X  I  [-  y  dV'  W£.(x"  -  Ax)</>2(x")  +  02]  -  I  dV"  Wl^ix"'  ~  Ax)<^,(x"')  +  02-]  I 

(3.11) 

Note  that  the  subscripts  1  and  2  on  the  wave  front  phase  (p  refer  to  the  source  of  the  wave 
front.  Performing  the  multiplication  indicated  in  equation  3.11  results  in 
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C,  = 

J  d^x  J  d^r  WIXS)WI{^' -  ^x)[4>x(S)<t>2{^')\ 

-j  d^x  Wl^{x)<j)i{x)a2 

-J  dV' Ax)</>2(^'>1 

+aia2 

+  I  dV  y  dV"  Wl^ix')W[^ix"'  -  Ax)  [Ux')Ux"')] 

-  J  dV  Wl^ix')M^)(^2- 

-  J  dV"  W[^{x"'  -  Ax)<l>2(x"')ay 

-\-oiiia2' 

-y  dV  J  d^x"  wi^ix')wi^(r  -  Ax)[MS')4>2{n] 

+  J  dV  Wl^{x)4>i{x)a2 

+  J  d^r  Wl^ix"  -  Ax)<i)2(x")ay 
—ai>a2 

d^S'''Wl^ix)Wl^ir'-Ax)[Mm2{S''')] 

+  J  d^xWl^{x)4>i{x)a2' 

+  J  d^x"'  Wl^{x"'  -  Ax)<t>2{x'")ai 

— Qia2'  (3.12) 

Equation  3.12  contains  the  random  quantities  <f)„  and  a„,  therefore  further  progress 
requires  a  statistical  approach.  Assuming  that  a„  is  a  zero  mean  random  process  indepen¬ 
dent  of  the  ensemble  average  of  equation  3.12  gives  the  desired  correlation  {C,} 

{C,)  =  J  d^x  J  d^x"  Wl^ix)Wl^ix"  -  Ax){MS)<t>2{n) 

-0 

-0 

+  (ol«2) 
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+  I  dV  y  Wl^{x')Wl^{r'  -  Ax){M^)US''')) 

-0 

-0 

+(o;i'0!2') 

-  y  dV  y  A-x''  -  Ax){U^')c}>,{x'')) 

+0 

+0 

-{avOi2) 

-  y  d^f  y  dV"VV»,(x)VF^^(x"'-Ax){<^i(f),/.2(x"')} 

+0 

+0 

-(ttia2')  (3.13) 


Wallner  has  shown  that  the  correlation  of  the  noise  process  a  is  given  by  (20) 

{anam)  =  fi{n  -  m)  (3.14) 

where  n  and  m  refer  to  the  sensors’  apertures,  is  the  variance  of  the  slope  measurement 
due  to  photon  noise,  and  6(x)  is  the  Dirac  delta  function.  Basically,  equation  3.14  states 
that  the  correlation  of  the  noise  processes  from  two  subapertures  is  zero  unless  the  sub- 
apertures  coincide.  Notice  that  none  of  the  noise  process  correlations  in  equation  3.13  are 
due  to  subapertures  that  coincide.  Therefore,  equation  3.13  reduces  to 

(C,)  =  y  d'*  y  dV'  Wl^{x)Wl^ix"  -  Ax){MS)US")) 

-^y  dV  y  A^x''' Wl^{x')W^x''' -  Ax){MS')US'")) 

-J  I  dV'VTl,(x')M^*.(x"-Af)(<A,(f')<A2(^")) 

-  y  dV  y  A^x'''  W[^{x)Wl^{x'''  -  Ax){Mx)U^''))  (3.15) 

Consider  the  third  and  fourth  terms  of  equation  3.15.  Their  forms  are  nearly  identical,  with 
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the  only  difference  being  the  fact  that  apertures  and  have  switched  positions.  It 
is  shown  in  section  3.3.3  that  the  correlation  of  the  phases,  {(pi{x)(f>2(x'))y  is  spatially  wdde 
sense  stationary,  with  the  value  of  the  correlation  depending  upon  only  the  magnitude  of 
the  aperture  separation.  Using  this  fact,  the  last  two  terms  have  the  same  value  and  may 
be  combined  so  (C,)  becomes 

{C,)  =  J  d^x  j  d-x''  Wl(x)Wl(r  -  ^x){4>,(x)US'')) 

+  I  dV  I  dV'VU»^(x')VU>^(f'"- 

-2 1  dV  J  d^r  W^x')Wl^(x''  -  Ax){(l>,{x')<i>2in)  (3.16) 

3.3.3  Evaluation  of  First  Term  in  (C,)  As  equation  3.16  shows,  the  first  term  of 
(C,)  is  given  by 

Term  I  =  J  d^x  J  d^x"  VUl,(5')VU£,(f"  -  Af)(</>i(x)</>2(f"))  (3.17) 

where  {4>i[x)<i>2{iP'))  is  the  correlation  of  the  two  wave  front  phases.  The  correlation  of  the 
phases  may  be  expressed  in  terms  of  the  phase  structure  function  as 

{0i(f)02(f"))  =  -\D,2ix,S")  +  i(<^?(f))  +  ^mx"))  (3.18) 

where  Z?i2(x,  it")  is  the  phase  structure  function  defined  by 

D,2(i,r)^{[M.x)-MS")f)  (3.19) 


Taking  note  of  the  fact  that  W’(x)  is  an  odd  function  and  that  <f>i  and  02  are  wide 
sense  stationary  random  processes,  the  integrations  which  result  when  equation  3.18  is 
substituted  into  equation  3.17  will  yield  zeroes  for  all  terms  in  equation  3.18  except  the 
term  containing  Di2ix,x").  So,  the  first  term  may  now  be  expressed  as 


Term  1  =  J  d^x"  Wl^{x)Wl^ix"  -  Ax)D,2{x,x")  (3.20) 


3-9 


Further  progress  requires  an  expression  for  Di2(x,x'').  For  spatially  stationary, 
isotropic  turbulence,  Lutomirski  and  Buser  have  shown  the  structure  function  for  the 
geometry  in  Figure  3.1  is  (12) 

Di2{x-x")  =  8TT^k^  j  dz'  J  dK\ 

x|l- Jo(A'x|(A^(^f)  +  (f-f")(|l-^)  )|  (3.21) 

where 

$„(A'x,A'„^')  is  the  psd  of  the  fluctuations  of  the  index  of  refraction,  n, 

A'x  is  the  wavenumber  component  perpendicular  to  the  z  axis, 

A'j  is  the  wavenumber  component  parallel  to  the  z  axis, 

Ap  is  the  vector  separation  of  the  sources, 

Jo  is  the  zero  order  Bessel  function, 

z'  is  the  distance  from  the  aperture  plane  to  some  point  along  the  2  axis, 

2,  is  the  altitude  of  the  point  sources, 

and  X  and  x"  are  points  in  the  aperture  plane  of  the  optical  system. 

Notice  the  psd  is  a  function  of  the  wavenumber  components  perpendicular  (A'x) 
and  parallel  (A'^)  to  the  2  axis.  By  making  the  assumption  that  is  separable  in  the 
perpendicular  and  parallel  wavenumber  components,  and  2',  Welsh  expresses  the  psd  as 
(22) 

<&„(A:x,  A',,2')  =  ^'„{Kj_,K,)Cl(z')  (3.22) 

Using  the  structure  function  given  by  Lutomirski  and  Buser  and  making  Welsh’s  assump¬ 
tion  of  psd  separability  allows  the  first  term  of  equation  3.16  to  be  written  as 

Term  1  =  J  d^x  J  dV' W’,(i)W;,(f"  -  Ai) 

xi-4ir^k^)  r  dz'  r  dA'xA'x$:.(A^x,0)C2(2') 

Jo  Jo 

X  |l  -  Jo  (a'x  |(A^  (0  -h  ix-x")  (1  -  0  (3.23) 
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Now  that  C\  has  been  introduced,  this  first  term  may  be  written  as 


Term  1  =  I  Az'  C‘^{z')wi{z')  (3.24) 

Jo 

where  Wi{z')  is  a  path  weighting  function  defined  by 

wi{z‘)  = -4Tr^P  j  d^f  j  &‘x"Wl^{x)Wl^{x"  -  Ax) 

dA-xA^x$UAx,0)|l- Jo(a’j.  (Ap)  | 

(3.25) 

The  goal  now  is  to  find  an  expression  for  tOi{z')  that  depends  upon  the  measurement 
geometry  alone  (i.e.  z,.  Ax,  Ap).  The  first  step  is  to  make  use  of  the  Kolmogorov  spectrum 
discussed  in  Chapter  2  to  assume  a  form  for  given  by  (9) 

KiKi,  A,)  =  0.033A'-^  (3.26) 

where  K  =  |Ax  +  A^J.  Substituting  equation  3.26  into  equation  3.25  yields 

wi{z')  = -4Tr^e  J  d^x  j  d^x"Wl^{x)Wl^{x"  -  Ax) 

X  j"  dh\  0.033A7"  |l  -  Jo  (Ax  [(A^  +  (x  -  x")  (^1  -  | 

(3.27) 

Now  consider  the  change  of  variables 

A'  =  A'x  |( A^  (0  +  (5  _  r )  -  0 1  (3.28) 

and 

dA"  =  (Ap)  +  (x  -  f")  dAx  (3.29) 
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Employing  the  integral  identity  (15) 


At  this  point,  the  path  weighting  function  depends  upon  the  measurement  geometry 
and  the  slope  sensor  aperture  weighting  functions.  Square  apertures  are  commonly  used 
in  wave  front  sensors,  and  the  transmittance  of  such  apertures  can  be  modeled  as 


^l{S)  =  ^rect  rect 


(3.32) 


where  rect  is  the  rectangle  function,  L  is  the  dimension  of  the  aperture  along  one  side, 
and  the  factor  ^  is  necessary  to  achieve  normalization  of  the  aperture  weighting  function. 
Assuming  that  the  slope  sensor  apertures  under  consideration  possess  a  transmittance 
function  of  the  form  given  in  equation  3.32  and,  without  loss  of  generality,  that  they 
measure  the  wave  front  phase  slope  in  the  x  direction  only,  then  the  gradient  of  the  slope 
sensor  aperture  weighting  functions  may  be  written  as 


wiAi) = if  K- + t)  -  *  ("  -  f )]  (i) 

where  6  is  the  Dirac  delta  function. 

Before  the  slope  sensor  gradient  functions  just  presented  can  be  useful  in  equation 
3.31,  the  absolute  value  quantity  in  equation  3.31  must  be  written  in  terms  of  the  x-  and 
y-directed  components  of  Ap,  x,  and  x".  The  absolute  value  term  can  be  seen  to  be 
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ApJ-)+{y-y")  fl--) 

\z,J  V  zj 


5 

6 


(3.34) 


where  Ap^  and  Apy  are  the  x-  and  2/-directed  components  of  Ap,  x  and  x"  are  the  x- 
directed  components  of  x  and  x",  and  y  and  y"  are  the  y-directed  components  of  x  and 
x".  It  is  convenient  to  represent  equation  3.34  with  the  function  f{x  —  x",  y  —  y") 

fix  -  x",y-y")  = 

Ap.  (I)  +  (X  -  X”)  (l  -  0 

(3.35) 


+ 


Ap,(0+(y-/)(l-0] 


Substituting  equations  3.35  and  3.33  into  equation  3.31  yields 


w,{z')  =  -lA&kH\^  j  d^f  j  d^x^'fix-x'^y-y") 

x{[0x+0-«(x-0],ect(0} 

X  {  [<  (x"  -  Ax  +  ^)  -  «  (x"  -  Ax  -  :^)]  reel  } 


(3.36) 


Taking  advantage  of  the  sifting  property  of  the  Dirac  delta  function,  the  integrations  over 
X  and  x”  yield 

Wi{z')  = -lA6k^Li‘*  J  dj/rect^^^  J  dy"rect^^  T 

X  {fi-Ax,y-y")  -  fi-Ax  -  L^y  -  y")  ~  f{-Ax  +  Li,y  -  y”)  A  f  {-Ax,y  -  y")} 

(3.37) 


At  this  point,  it  is  convenient  to  define  the  function  g{y  -  y")  as  the  bracketed  term  of 
equation  3.37 
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9{y-y")  =  f(-Ax,y-  y")  -  f{-Ax  -  Li,y-  y") 
-/(-Ax  +  Li,y-  y")  +  f  {-Ax,  y  -  y") 


(3.38) 


Use  of  equation  3.38  allows  the  path  weighting  function  to  be  written  as 

wi{z')  = -lAGk^Li'^  J  dr/rect  j  dy"  reel  ^  g{y  -  y" 


)  (3.39) 


Making  the  change  of  variables  u  =  y  —  y"  and  v  =  allows  equation  3.39  to  be  written 


Wi{z')  =  —lA6k^Li‘^  j  durect 


^)/ 


Z' V  —  ~  ~  Aw\ 

durect  ^ ^ - j  g(u)  (3.40) 


One  more  change  of  variables,  letting  u'  =  u  +  f ,  results  in 

Wi{z’)  = -lAdk^Li"^  j  dug{u)  j  du' rect  j  ^ 


v'  -  u  -  Ay 
L'l 


Notice  that  the  integral  over  v'  is  simply  the  autocorrelation  of  a  rect  function.  Thus 
evaluating  the  integral  over  v'  yields 


Wi{z')  =  -lA6k^  j  dMg(M)  tri  ^ 


(3.42) 


where  tri  is  the  triangle  function  defined  by 


tri(x)  = 


1  -  |x|  -1  <  X  <  1 

0  elsewhere 


(3.43) 


Evaluating  g{u)  through  use  of  equations  3.38  and  3.35  and  placing  the  result  in  equation 
3.42  yields  the  expanded  expression 

Wi{z')  =  -lA6k^L'[^  J  du  tri 

{ h'  (i)  '('■?:)  '^"1  h'  (0 + "  -  J;)] } 
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pression  is 


The  weighting  function  in  equation  3.47  above,  coupled  with  equation  3.24,  completes 
the  evaluation  of  the  first  term  of  {C,).  The  second  term  in  equation  3.16  will  be  examined 
next. 


3.3.4  Evaluation  of  Second  Term  in  {C,)  The  second  term  of  equation  3.15  is  given 

by 

Term  2  =  J  d'^x'  J  d'f"'  -  Ax){<j>,{x')MS"'))  (3.48) 

Notice  that  this  term  is  identical  in  form  to  the  first  term  with  one  exception:  Li  has  been 
replaced  by  £3  in  the  aperture  weighting  function.  Performing  the  same  analysis  for  this 
term  as  that  performed  for  Term  1  will  yield  a  result  identical  to  that  obtained  for  Term 
1,  but  with  Li  replaced  by  L2.  Thus,  the  second  term  can  be  written 


Term  2  =  /  d^'  C^(z')  ^2(2') 

Jo 


where  W2(z')  is  given  by 


(3.49) 
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Now  that  expressions  for  the  first  two  terms  in  equation  3.16  have  been  obtained,  it 
is  time  to  evaluatf*  the  third  term. 

3.3.5  Evaluation  of  Third  Term  in  (C^)  The  third  term  in  equation  3.16  is 

Term  ^  J  J  dV'  (3.51) 

This  term  is  essentially  the  cross  correlation  of  the  phase  slope  measurement  on  an  aperture 
of  dimension  Li  with  the  phase  slope  measurement  on  an  aperture  of  dimension  L2,  where 
the  two  apertures  are  separated  by  Ax.  The  derivation  of  an  expression  for  the  third 
term  follows  the  same  development  as  that  for  the  first  term  until  models  are  assumed 
for  the  gradients  of  the  slope  sensor  weighting  functions.  Thus,  the  starting  point  of  this 
derivation  is 

Term  3  =  dz'  C^(z')  Wsiz’)  (3.52) 

Jo 

where  1^3(2')  is  given  by 

W3{z')  = -2.92k^  J  dV  J  d^x"Wl^{iT)Wl^{x"  ~  Ax) 

X  |(A^  +(x'-  x")  (1  -  7)  '  (3-53) 

Once  again,  the  slope  sensors  are  assumed  to  measure  the  wave  front  slope  in  the  a:-direction 
only,  so 

=  ;^  [<(“’'  +  t)  - '  (*'  - 1)]  it) 

where  n  =  1  or  2.  The  absolute  value  quantity  present  in  equation  3.53  is  again  designated 
with  the  function  /  defined  in  equation  3.35.  Substituting  equations  3.35  and  3.54  into 
equation  3.53  and  performing  the  integrations  over  x'  and  x"  yields 


3-17 


(3.55) 


W3{z')  =  —2.92k^L^^L2^  J  dy'rect^^^  J  dy"  reel 

X  {/ (-Ax  -  -!,")-/ (-Ax  -  j”) 

-f  (^-Ax  +  :^  +  |-,  j  -  j”)  +  /  (-Ax  +h_h^y_  j,»)| 

As  was  done  during  the  derivation  of  Term  1,  let  the  bracketed  term  be  designated 


93iy'-y")  = 

/  (-Ax  -  -  j,")  -  /  (-Ax  -  -  |.,9  -  j,”) 

-/  (-Aa:  +  y  +  Y'^~  ^  ^  +  T  "  Y’ ^ 

Use  of  equation  3.56  allows  equation  3.55  to  be  written  as 

Wsiz')  =  -2.92k^L';'^L2^  J  dy'rect(^|-^  J  dy"  rect  gsiy' -  y”)  (3.57) 

Making  the  same  substitutions  of  variables  as  accomplished  in  subsection  3.3.3  results  in 

waiz')  =  ~2.92k^L'[^L2^  j  du  g3{u)  J  dn'rect  ^ ^ — —'j  (3.58) 

The  integral  over  v'  in  equation  3.58  is  simply  the  correlation  of  two  different  sized  rect 
functions.  Performing  the  correlation  yields 

/  du'rect  (£)  rect  = 

u  +  +  Ay  for  -  (^^^  -f-  Ay)  <u<  “  Ay 

Ly  for  <  u  <  kz.-i'i  _  (3.59) 

-  u  —  Ay  for  ^a=ia.  _  <  u  <  -  Ay 

0  otherwise 

Once  again  assuming  Ay  =  Apy  =  0,  use  of  equation  3.59  in  equation  3.58  results  in 
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W3{z')  = 


-2.92k^L^^L2^  {  I  dti  (tt  + 


^2  +  ^1 


lfif3(u)  +  ii  duff3(u) 


+ 


LL  ( 


Z/j  +  ^2 


u  g^iu) 


(3.60) 


where  gsiu)  has  evolved  to  the  following  expression 

Z>2  —  L\ 


fl'3(w)  = 


Ap: 


-{[ap,|+(-Ax-  ^ 


2 

Li  +  Z/2 


)O-0] 


+ 


A  (a  -^1  -^2 

Ap;r—  +  (  -Ax  +  - - - 

2,  \  2  . 


-{h0+(- 

+  {  [Ap.i  +  (-Ax  +  (l  -  0]  + 


(3.61) 


Letting  7  =  Apx  -  Ax  and  performing  some  algebraic  expansions,  gsiu)  can  be  written 


as 


53(w) 


■( 


2  2,  V  2 


(ax -i^)r  001-0] 


2'>  6 


-I 

-{ 


L1  +  L2  z'  f ^  Zj  +  ZzAT  {  ^'Afl 

Li+L2,z'f  Zi  +  Z2A1^,  f  f,  z''' 

+  +K'-5:. 


7  + 


Li  —  L2  ,  z' 


+  £(Ax-ii0:j)]%K‘-0]y  (“2) 


The  first  step  in  making  the  integrals  over  u  dimensionless  is  to  let  R  be  the  ratio  of 
The  result  is 

L'i 


103(2')  = 


{ (“  ~T~~)  A-  -' 


2.92k^R-^L-'' 


3 

r  Lj(B-i) 
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I>2(W+U  , 

I  ^  ( 


L2iR+l) 


-  W  53(w) 


(3.63) 


where 


,  ,  f  r  L2(1-  R)  Z' 

*(“)=  ’  + - 2 - +r 


L2{R+1)  I  2' 

(ax '))]’+[„( 

i-i 

2  2, 

V  2  yj  L  V 

z. 

L2{R+1)  z' 

z' 

1  — 

2  z. 

V  2  ;J  [  V 

L2{R-1)  z' 

(Ax-"^(V”)l“+k 

1-i: 

2  z. 

V  2  yj  L  \ 

2, 

Making  the  change  of  variables  u'  —  reduces  ^3(2')  to 
^3(2')  = 


(3.64) 


(3.65) 


where 


gsiu') 


'  \  1  21  6 


-{I5 

•Hi 


J-  +  (^+1)  +^(^_  (i?+  1) 
L2  2  z,  \  L2  2 


)]’^K-0]T 


(3.66) 


This  completes  the  evaluation  of  the  third  term  in  equation  3.16.  The  third  term  is 
given  by 

Term  3  =  J  dz' Cl{z')w3{z')  (3.67) 
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where  ^3(2')  is  given  in  equation  3.65  above.  The  three  terms  in  equation  3.16  are  combined 
in  the  next  section. 


8.3.6  Combination  of  Terms  Equation  3.16  has  been  broken  into  three  terms  and 
each  term  evaluated  separately.  The  three  terms  may  now  be  combined  to  yield 


{C,)  =  I  dz'CfXz')wc{z') 


(3.68) 


where  wc{z')  =  Wy{z')  +  W2{z')  —  w^^z').  Using  the  results  derived  previously,  and  letting 
R=  Wc{z')  may  be  written  as 


I  yl 

Wc{z')  = R~^  J  du' tri(«') 


'Hte-oa’-K'-si 

-{ 


/  \  1  2  'i  6 


'RL2~^^  J 


-1.46k^L2^  J  du' tri(u') 


2 

1  ,  z'V 

'1 

+ 

u'il - 

> 

.  V  zj. 

J 

l2 

+ 

-( 


/  \  1  2  1  6 


-\-2.92k^R-ni 


K'-0  , 

V  2 


IL2  z,  \  L2 


7  .  -2 
7-  +  1  + 

L2  z. 


0-«) 


du'  ^u'  +  ^  ff3(u')  +  du'  flf3(u') 


(3.69) 


where 
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[[7  (1-^)  z'  (^x  0  -  R)\Y  \,f  z' 


\-L  _  (-^  +  ,  (^+  1) 

.Z2  2  z,  \  L2  2 

'JL  +  (^  +  1)  I  _  (^+1) 

.Z2  2  2,  \  Z(2  2 


+  [1;+ 


7  (i?  -  1)  2'  /Ax  (i?-  1) 


U2 


Equation  3.69  is  the  most  general  form  of  wc{z').  However,  the  path  weighting  func¬ 
tion  can  be  simplified  further  if  the  assumption  is  made  that  1  —  ^  ss  1.  This  assumption 
is  often  valid  when  the  point  sources  are  artificial  guide  stars  created  in  the  mesospheric 
sodium  layer.  In  this  case  z,  is  on  the  order  of  90  km,  and  z'  is  often  no  greater  than  20 
km  because  «  0  for  z'  >  20km.  The  simplified  form  of  the  path  weighting  function  is 

wc{z')  = R~^  j  du' tri(u') 

X  [(^)’ + “i  ‘-[{-k- ^  “i '  -  [(sk + ')’ + “i ' ) 

—  1.46A:^T2’  /  du'tri(u') 

’‘H(i)  ■^“i  +“'i  +"i  1 


2  1  6 


+  1  +  «'■ 


2  1  6 


+2.92k^R-H-^  [/-(^)  r  + 

,/^/(i?+l) 


where 


= { [i + + '“'if  -  { [i  -  + f 


(3.72) 
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The  simplified  path  weighting  function  given  by  equation  3.71  is  examined  in  the  next 
section.  Recall  that  the  z'  dependence  of  Wc{z:'}  is  contained  in  7,  where  7  :=  —  Ax 

and  may  be  interpreted  as  the  distance  between  the  optical  beams  as  a  function  of  altitude 
(see  Figure  3.1).  In  the  lesults  which  follow,  Wc(x')  is  plotted  versus  7  rather  than  2'. 

3.4  Form  of  C\  Path  Weighting  Function 

The  form  of  the  path  weighting  function  for  the  new  measurement  technique 
is  one  of  the  primary  objectives  of  this  research.  As  already  shown  in  equation  3.68, 
there  is  an  integral  relationship  between  the  measured  wave  front  phase  correlations  and 
C^.  If  the  path  weighting  function  Wc(z')  approximates  a  Dirac  delta  function,  then  the 
measured  quantity  (C,)  can  be  directly  related  to  A  normalized  plot  of  the  simplified 
path  weighting  function  given  by  equation  3.71  is  shown  in  Figure  3.3  for  three  reasonable 
values  of  the  ratio  R  = 


0  2  4  6  8  10 


Y/L^ 

Figure  3.3.  Normalized  Path  Weighting  Function 


The  weighting  function  displays  the  desired  characteristic  of  being  sharply  peaked  at 
the  altitude  where  the  optical  beams  intersect  (7  =  0)  and  decaying  rapidly  to  zero  away 
from  the  intersection  altitude.  The  plot  is  shown  one  sided  because  the  path  weighting 
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function  is  symmetrical  about  the  intersection  altitude.  Notice  that  the  7  values  on  the  x 
axis  of  Figure  3.3  are  normalized  by  h.  As  an  example,  if  /,[  =  10  cm  and  R  -  0.9,  then 
the  path  weighting  function  is  zero  by  the  time  7  =  20  cm.  The  figure  also  demonstrates 
the  role  of  iZ  =  ^  in  shaping  the  path  weighting  function.  j\s  R  I,  the  path  weighting 
function  becomes  more  sharply  peaked.  This  makes  sense  because  as  R  approaches  one, 
higher  frequencies  are  being  removed  from  the  sensor  measurements.  In  other  words,  the 
pass  band  of  the  filtering  operation  accomplished  by  the  subtraction  of  the  measurements 
is  becoming  more  narrow.  This  corresponds  to  reducing  the  eddy  sizes  which  contribute 
to  the  correlation.  Such  eddies  simultaneously  intersect  both  beams  only  when  7,  which 
is  the  separation  of  the  two  beams,  is  small.  Therefore,  as  the  eddy  sizes  contributing  to 
the  correlation  decrease,  the  correlation  decay  rate  increases. 


3.5  Resolution  of  Path  Weighting  Function,  Wc{z') 

As  Figure  3.3  shows,  the  path  weighting  function  is  not  a  Dirac  delta  function,  but 
does  have  its  peak  value  where  the  beams  intersect  and  decays  to  zero  away  from  the 
intersection  point.  The  width  of  Wq  with  respect  to  z  may  be  considered  a  measure  of  the 
vertical  resolution.  The  narrower  Wc{z)  is,  the  better  its  vertical  resolution  will  be,  and 
thus  a  more  accurate  estimate  of  C^iz)  will  be  obtained. 


A  method  similar  to  that  presented  by  Rejack  (17)  is  used  to  measure  the  vertical 
resolution.  Recall  from  Figure  3.1  that  the  transverse  distance  between  the  two  optical 
beams  is  given  by 

z' 

7  =  Apj, - Ax  (3.73) 


where  Ap^.  is  the  source  separation,  z,  is  the  source  altitude,  and  Ax  is  the  sensor  separa¬ 
tion.  The  altitude  Zq  where  the  optical  beams  cross  is  found  by  setting  7  =  0  and  solving 


for  Zo 


Ax  z. 


(3.74) 


Now  let  the  measure  of  the  width  of  Wc(z')  be  the  point  7  =  7'  where  7'  is  the  one 
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sided  root  mean  square  (rms)  width  of  Wc  defined  by 

/o”  “’c  (fr) 

Denote  the  altitude  corresponding  to  7'  as  where 

(7'  +  Ax)2, 


(3.75) 


Apx 


(3.76) 


Due  to  the  even  symmetry  of  the  path  weighting  function  about  the  intersection  altitude, 
the  width  of  Wc  with  respect  to  z  is 


AZ  =  2(2:0  —  2y  )  = 


2Y  Z. 


(3.77) 


Let  0  denote  the  angle  between  the  optical  beams.  Referring  to  Figure  3.1,  basic 
geometry  yields 

(jt:^)  (3-78) 


0  (  Apx 

-  =  arctan  I  — - - 

2  V2(z,  -  2:0) 


The  assumption  that  z,  >  2o  allows  0  to  be  equated  to  the  angular  separation  of  the 
sources  given  by 


0 

-  =  arctan 
2 


(t) 


Solving  for  Apx  yields 


Apx  =  22:,  tan 

Substituting  equation  3.80  into  equation  3.77  results  in 

-v' 

A2:  = 


(3.79) 


(3.80) 


tan  (|) 


(3.81) 


All  the  expressions  necessary  to  calculate  the  vertical  resolution  of  wdz')  are  now 
present.  An  example  calculation  can  be  performed  based  upon  the  results  in  Figure  3.3. 
For  R  =  0.9  and  Li  =  10  cm,  the  rms  width  of  wc{zo)  is  7'  =  0.136  m.  Using  this  value 
in  equation  3.81  yields  the  vertical  resolution  as  a  function  of  0  plotted  in  Figure  3.4. 
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9  (degrees) 


Figure  3.4.  Vertical  Resolution  of  wc 


The  figure  demonstrates  the  relationship  of  the  vertical  resolution  to  the  separation  of 
the  point  sources.  As  the  point  source  separation  increa.ses,  so  does  the  vertical  resolution. 
The  reason  for  this  behavior  is  that  a  larger  source  separation  increases  the  angle  at  which 
the  beams  intersect.  As  a  result,  less  vertical  movement  away  from  the  intersection  altitude 
is  needed  to  reach  a  point  where  the  separation  of  the  beams  is  greater  than  the  size  of 
the  eddies  which  are  contributing  to  the  measurement  correlations. 


3.6  Derivation  of  V  Path  Weighting  Function 

3.6.1  Overview  The  second  objective  of  this  research  is  to  extend  the  new  mea¬ 
surement  technique  to  the  remote  sensing  of  transverse  atmospheric  winds.  Accordingly, 
the  goal  of  this  derivation  is  to  relate  the  derivative  of  (C,(r))  with  respect  to  r  to  the 
vertical  wind  profile  V  through 


^  ^([^i(0  -  ar(0][g2(<  +  r)  -  sr(t  +  r)]) 


I  dz'Cl{z')V{z')-wviz') 


(3.82) 

where  (C'(r))  is  the  partial  derivative  with  respect  to  r  of  {C,{t)),  z'  is  a  position  along 
the  optical  path,  and  Wv{z')  is  the  wind  path  weighting  function.  This  derivation  closely 
parallels  the  Cl  path  weighting  function  derivation  just  completed  and  will  often  refer 
to  it.  Eventually,  the  path  weighting  function  is  reduced  to  an  expression  that  depends 
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only  upon  the  mecisurement  geometry.  The  reasons  for  taking  the  partial  derivative  with 
respect  to  r  and  evaluating  it  at  r  =  0  will  be  discussed  at  the  appropriate  points  within 
the  derivation.  The  first  step  in  this  derivation  is  to  find  an  expression  for 

3.6.2  Slope  Signal  Time-lagged  Correlation  As  shown  above,  the  single  realization 
time-lagged  correlation  is  defined  as 

C,{t)  =  [6i(<)  -  5i.(t))[s2(<  +  r)  -  S2'it  +  r)]  (3.83) 

where  s„{t)  is  the  output  of  slope  sensor  n  at  time  t  and  r  is  the  time  delay  between 
measurements.  Using  the  wave  front  slope  sensor  output  model  given  by  equation  3.10, 
equation  3.83  can  be  written  as 


C,{T)  = 
{[-/ 


I  [-  I  d  V'  W[^(x"  -  ^  +  r)  +  a,(t  +  r)] 


-  -J  dV"  Wl,{r’  -  Ax)Mx"',  /  +  r)  +  arit  +  r)  | 


(3.84) 


where  a„(t)  is  the  slope  measurement  noise  from  sensor  n  at  time  t  and  (j)„{x,t)  is  the 
wave  front  phase  arriving  from  source  n  at  time  t.  Performing  the  multiplication  results  in 


C,{t)  =  I  d^x  J  d^x''Wl^{x)Wl^{x''-Ax)[MS,t)Mi'',t  +  T)] 

-J  d'^^Wl^{x)(l)i{x,t)a2{t-\-T) 

-  j  dV'  Wl^{x"  -  Ax)<^2(x",<  +  r)a,(0 
+ai(t)Q2(t  +  t) 

+  |dV|  dV"VT’^(f')lT’,(f"'-Ax)[<^i(f',0<^2(^",f  +  r)] 

-  J  d^x^Wl^ix')MS',t)aAt-\-T) 

-  J  d^x'"  -  Ax)<i>2{S''\t-\-T)aAt) 

+ai>(<)a2'(<  +  7") 
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-y  dV  y  d-x''  Wl^{x')Wl^(x''  -  Ax)[M^,t)(l>n(x'\t  +  T)] 

+  j  d^x"  Wl^{x)(l)i{x,t)a2{t  +  t) 

+  j  d'^x"  Wl^{x"  -  Ax)<j>2(x'',t  +  T)ai<{l) 

-av{t)a2(t  +  r) 

-J  d-x  I  d^r'w[^ix)wi^(r'-Ax)[Mx,t)M^'\i  +  T)] 

+  J  d'^xWl^{x)(f)i{x,t)a2-it  +  T) 

+  J  d^f"'  W^x"'  -  Ax)<l>2{x"',t  +  T)a,it) 

— ai(i)a2<(<  +  r)  (3.85) 

Once  again,  further  progress  requires  a  statistical  approach.  Assuming  that  a„  is  a 
zero  mean  random  process  independent  of  4>n,  and  using  the  fact  that  the  noise  processes 
from  separate  apertures  are  uncorrelated  (see  equation  3.14),  the  ensemble  average  of 
equation  3.85  is 

(C.(r))  =  J  d^x  J  d^x''  Wl^{x)Wl^ix''  -  Ax){MS,t)M^\t  +  T)) 

+  Jd^S'j  dV"Vy2,(f')M^/,(f'"-Ax)(</>i(x',<)<^2(f'",<  +  r)) 

-y  dVy  dV'V^£,(f')W^^,(f"-Ax)((^,(f',t)<^2(r,t  +  r)) 

-J  d^f  y  dV'W^/,(f)W^‘,(x'"-Af)(<^i(x,f)</>2(f"',t  +  r))  (3.86) 

The  third  and  fourth  terms  of  equation  3.86  can  be  combined  for  the  same  reasons  that 
the  third  and  fourth  terms  of  equation  3.15  were  combined.  Thus,  the  ensemble  average  is 
given  by 

(C,(r))  =  y  d^x  y  dVvy^,(x)iy'.(x"- Ax)(<^i(x,0^2(?",<  +  r)) 

+  y  d^x'  y  dV"iy^^(x')iy/,(x"'-Ax)(<^i(x',<)<A2(f'",f  +  r)) 

-2 y  dV  y  d^x^'Wl^ix')Wl^{r-Ax){Mx',t)US'\t  +  T))  (3.87) 

Equation  3.87  consists  of  three  terms.  Each  of  these  terms  is  evaluated  separately 
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and  the  results  combined. 


3.6.3  Evaluation  of  First  Term  in  {C,(r))  The  first  term  of  equation  :5.87  is  given 

by 

Term  1  =  J  d'x  j  dV' -  Af  <  +  r ))  (3.88) 

where  {(j)i(x,t)(f)2{x" ,t  +  r))  is  the  cross  correlation  of  the  wave  front  phases.  Once  again, 
the  cross  correlation  of  the  wave  front  phases  corresponds  to  the  structure  function  as 

(</),  (f,  t  )d)2(f",  <  +  r)}  -  -~D,2ix,  x'\t)  +  t))  +  t  +  r))  (3.89) 

where  the  structure  function  is  now  defined  as 

Du{x,  f",  r)  =  ([<^i(f,  0  -  d>,(f",  t  +  r)]^)  (3.90) 

Due  to  the  odd  symmetry  of  IT’(f)  and  the  fact  that  <^i  and  4>2  are  wide  sense  stationary 
processes,  only  the  term  containing  Di2{x,x",t)  will  yield  a  non-zero  result  when  equation 
3.89  is  placed  in  equation  3.88.  An  expression  for  Di2{x,x")  is  given  in  equation  3.21. 
However,  equation  3.89  requires  that  the  time  lag  be  accounted  for  in  the  structure  function. 
This  is  accomplished  in  the  same  manner  used  by  Lee  and  Harp,  which  is  to  replace 
(f  -  x")  in  equation  3.21  with  (x  —  x")  ^1  -  -  V{z')t.  This  substitution 

accounts  for  the  distance  traveled  by  the  turbulent  eddies  between  measurements.  Figure 
3.5  illustrates  the  situation.  Recall  that  the  measurement  of  (f>2  is  made  r  seconds  after 
the  measurement  of  </>i.  The  eddy  present  at  a  point  z'  along  the  ray  path  from  p2  at  time 
t  was  a  distance  -V{z')t  away  at  time  t  -  t.  Thus,  the  time-lagged  measurements  are 
equivalent  to  making  simultaneous  measurements  for  the  geometry  depicted  with  dotted 
lines  in  Figure  3.5.  Making  such  a  statement  requires  use  of  Taylor’s  frozen  flow  hypothesis, 
which  assumes  that  the  eddies  do  not  decay  as  they  are  blown  across  the  optical  line  of 
sight  (19).  So.  making  the  substitution  into  the  structure  function  given  by  Lutomirski 
and  Buser  yields 
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Altitude  z 


Figure  3.5.  Effect  of  V  on  Measurement  Geometry 


Di2(x  -  x",r)  =  Stt 


X 


J  dz'  I  dKj^MK\,0,z') 


(3.91) 


where,  as  before, 


^n(^jL>  K^,z')  is  the  psd  of  the  fluctuations  of  the  index  of  refraction, 
Ap  is  the  vector  separation  of  the  sources, 

Jo  is  the  zero  order  Bessel  function, 

z'  is  the  distance  from  the  aperture  plane  to  some  point  along  the  z  axis, 
Zf  is  the  altitude  of  the  point  sources, 
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and  X  and  x”  are  points  in  the  aperture  plane  of  the  optical  system. 

Substituting  equation  3.91  into  equation  3.88  and  using  equations  3.22  and  3.26  to  describe 
the  psd  of  the  fluctuations  of  the  index  of  refraction,  the  first  term  in  equation  3.87  is 


Term  1  = 


dz'C^iz')  j  d"£  J  d^x"  Wl^ix)Wl^ix"  -  Ax) 


dK,  0.033/t'; 


X  n  -  Jo  A' 


{Ap)  (^0  +{x-  x")  (l  -  0  -  y{^')r  )  I  (3.92) 


Performing  the  integration  over  K±  in  a  similar  fashion  to  the  integration  shown  in  equa¬ 
tions  3.28  through  3.30  yields 


Term  1  = -1.46P0  dz'C'Cz')  J  d^x  J  d'f"  VP£^(f)VT’^(x"  -  Af) 

(A,J)(0  +  (*-i") (l-0 -{’(*'> 


(3.93) 


which  may  also  be  written  as 


Term  1  =  -lA6k^  dz' Cl{z')  j  d^x  J  dV' lT0(f)fT0(f"  -  Ax) 

5 

X  I  +  (a:  -  x")  -  0  -  14(^')7’]  +  ^Py~  +  {y  -  y”)  (^1  -  0  -  v;,(2')7-  I 


(3.94) 


where  and  Vy  are  the  x-  and  7/-directed  components  of  V,  Apx  and  Apy  are  the  x- 
and  j/-directed  components  of  Ap,  x  and  y  are  the  x  and  y  directed  components  of  x,  and 
x"  and  y"  are  the  x  and  y  directed  components  of  x".  To  facilitate  a  later  step  where 
the  above  equation  is  differentiated  with  respect  to  r,  some  algebraic  manipulation  of  the 
squared  terms  results  in 
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Term  1  = -1.46P  dz' C^z')  j  d^x  j  dV' -  Af) 

X  I  +  {x  -  x")  ^1  -  +  V^(z')t^  -  2V^{z')t  ^Pxy  +  (a:  -  x")  ^1  - 


^Py—  +  iy-  y")  1  -  — 


+  V;-{z’)t^  -  2V^{z')r 


^Py—  +  {y 


(3.95) 


For  now,  Term  1  will  be  left  in  the  form  given  by  equation  3.95.  Differentiation  of  this  ex¬ 
pression  with  respect  to  r  is  easily  accomplished,  so  further  simplification  is  not  necessary. 


3.6.4  Evaluation  of  Second  and  Third  Terms  in  {C,{t))  Notice  that  the  second 
and  third  terms  in  equation  3.87  have  the  same  form  as  Term  1  with  two  exceptions: 
they  contain  different  dummy  variables  of  integration,  and  they  contain  different  aperture 
weighting  functions.  The  different  dummy  variables  of  integration  are  of  no  consequence. 
Additionally,  no  information  about  the  aperture  weighting  functions  in  Term  1  was  nec¬ 
essary  to  obtain  equation  3.95.  Thus,  the  second  and  third  terms  can  both  be  expressed 
in  a  form  identical  to  equation  3.95,  but  with  replaced  with  the  appropriate  aperture 
weighting  function  gradient.  The  second  term  in  equation  3.87  is 


Term2= dz' C^z')  J  dV  J  d^x"' Wl^ix')Wl^{x"' -  Ax) 

X  I  +  ix'  -  x'")  (|l  -  f )]'  +  '  2V,iz’)T  ^Aprj  +  (x'  -  x'")  (|l  - 

-1  [ap,-  -1  iy'  -  n  (1  -  -)1'  +  V^Hz')r^  -  2V^{z')t  Lpy-  +  {y'  -  y"')  (l  - 

\  /  J  L  \  / 


(3.96) 


and  the  third  term  is 

TeTm3  =  2.92k^  r  dz'Cliz')  f  d'x'  /  dV' VF’^(x')VT'_(x"  -  Ax) 


0 

Ap.-  +  (x'  -  x")  fl  -  -h  Vf{z')T^  -  2VAz')t 

Apr—  +  (x'  -  x")  ( 

‘-01 

1 

1  z,  \  Z,Ji 

L  2,  V 

Z,J\ 
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+ 


Ap,-  +  {y'  -  y")  fl  -  +  Vyi^'W  -  ‘^Vy{z')r  [Ap„f  +  {y'  -  y")  (l  -  f ) 

\  ■2’s  /  J  L  \  / . 


(3.97) 


These  three  terms  will  be  combined  and  the  result  differentiated  with  respect  to  r  in  the 
next  section. 

3.6.5  Differentiation  of  {C,(r))  With  Respect  to  t  The  previous  section  has  pro¬ 
vided  equations  for  the  three  terms  that  compose  (C,(t)}.  The  next  step  in  arriving  at  an 
expression  with  the  form  given  by  equation  3.82  is  to  combine  the  three  terms  and  differ¬ 
entiate  the  sum  with  respect  to  r.  Setting  r  =  0  in  the  result  of  the  differentiation  yields 
an  expression  of  the  desired  form.  Before  taking  these  steps,  however,  some  justification 
of  these  steps  is  necessary. 

One  of  the  first  techniques  for  sensing  V  using  crossed  optical  beams  was  to  measure 
the  difference  in  the  time  of  arrival  of  a  turbulent  eddy  at  one  beam  and  its  time  of  arrival 
at  a  second  beam.  This  time  delay  r  can  be  measured  by  making  time-lagged  correlations 
of  the  optical  beams’  wave  front  phases.  The  correlation  is  calculated  for  several  values  of  r, 
and  the  value  of  r  which  yields  the  peak  correlation  value  is  used  to  compute  V .  However, 
Lawrence  et  al  have  pointed  out  that  this  method  gives  results  which  are  systematically 
high,  particularly  for  low  wind  velocities  (10).  At  low  wind  velocities,  eddy  decay  causes 
the  turbulence  crossing  the  second  beam  to  bear  little  resemblance  to  that  crossing  the 
first  beam.  Therefore,  the  correlation  peak  occurs  sooner  than  it  should,  and  an  alternate 
approach  is  needed. 

This  problem  can  be  avoided  by  measuring  the  slope  of  the  correlation  function  at 
r  =  0.  In  practice,  the  slope  can  be  measured  by  subtracting  the  correlation  value  at  time  t 
from  the  correlation  value  at  time  t  +  r  and  dividing  the  result  by  r.  The  time  difference  r 
can  be  sufficiently  small  that  there  is  little  eddy  decay  between  measurements.  Evaluating 
the  slope  at  r  =  0  results  in  an  expression  with  the  form  of  equation  3.82.  Analytically, 
differentiating  the  time-lagged  correlation  function  with  respect  to  r  yields  an  expression 
for  the  slope,  and  setting  t  —  0  evaluates  the  slope  at  zero  time  lag. 
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Proceeding  with  the  derivation,  terms  1,  2,  and  3  are  combined  to  yield 


(C,(Af,r))  = 


-lA6k^  dz'Cliz')  J  d"f  y  Af) 

X  I  +  (x  -  x")  +  V^(z')t^  -  2V^(z')r  +  ix-  x")  ^1  - 

+  ^Py~ (y  ~  y"'^  ~  +V^(z')t^  -2Vyiz')T  Apy^  +  {y-y")^l- I 

-1A6P  dz'Cliz')  y  d'*x'  y  d^x"' VP*^(x')VP’^(x"' -  Af) 

X  I  +  (x'  -  x'")  (^1  -  f )]'  +  Vx(^'y  -  2V;(^')^  ^Prj  +  {x'  -  x'")  (^1  -  pj 

+  [ap,^  +  {y'  ~  y'")  (l  -  0]'  +  V^\z'y  -  2Vy{z')T  ^Ap^  +  (y'  -  y'")  (l  -  0]  | 

-\-2.92P  p'  dz' Cl{z')  y  dV  J  dV' VF£,(x')H^£,(x"  -  Ax) 

X  I  [ap,^  +  (x'  -  x")  (^1  -  0]'  +  VxVy  -  2K(2')^  [^pX  +  ix'  -  x”)  (l  -  0 

+  Apy^  +  (y'  -  y")  f )  +  V^\z'y  -  2Vy{z')T  Apy-  +  (y'  -  y")  (l  -  -)]  | 

\  /  J  I  \  ^5  /  J  1 


(3.98) 


The  partial  derivative  with  respect  to  r  of  equation  3.98  is 

a(C,(Ax,r))  ^ 
dr 

-IA6P  p’  dz'Cliz')  J  d^x  j  dV'tV£,(x)VP£,(x"- Ax) 

x^  {  +  {x  -  x")  ^1  -  0j  +  V^{z')t^  -  2V^{z')t  Api^  +  (x  -  x")  ^1  -  0 

+  [Apy^  +  (y-y")(l-0]  +V^(z')t^  -2Vy{z')T  Apy^  +  (y-y")^l-0  | 

X  ^2V^\z'}r  -  2Kiz')  [ap0  +  (x  -  x")  (^1  -  0 

+2P/(^')r  -  2Vy{z')  [ap,^  +  (y  -  y")  (l  -  0]  J 
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5 

^6 


+ 


X 


-lA6k^  dz' Cl{z')  J  dV  J  dV"  -  Af) 

I  +  {x'  -  x'")  ^1  -  +  V;{z')t^  -  2V^{z')t  +  (a;'  -  a;'")  ^1  - 

Ap,-  +  {y'  ~  y'")  fl  -  f )]'  +  V;^(^')r^  -  2Vy(z')T  [ap,-  +  (p'  -  y'")  (l  -  • 
\2V^{z')t  -  2Kiz')  [ap.  +  (a:'  -  x"')  (|l  - 
+2V^Hz')t  -  2V,{z')  [ap,^  +  (y'  -  y'")  (l  -  0]  } 

+2.92k^  dz'Cliz')  J  d"f'  J  d^x^'  W[^{x')Wl^(x"  -  Ax) 

X^  I  [ap^^  +  (ar'  -  a;")  ^1  -  0j  +  V^{z')t'^  -  2V^{z')t  Ap^^  +  (x'  -  a;")  ^1 

+  [ap,-  +  {y'  -  y")  (l  -  -)!'  +  V^{z')T-^  -  2V,(z')r  fAp,f  +  (p'  -  y")  fl  -  - 
X  |2Vj(2')r  -  2V.(J')  [ap,^  +  (i'  -  x")  (l  - 
+2V/(x')r  -  21^, (z')  [Ap,i  +  (p'  -  y")  (l  -  0]  } 


(3.99) 


Setting  r  =  0  gives 


aC,(Af,r), 
dr 


T  =  0  - 


-\Ade  p  dz'C^(z')  J  d^f  j  d^x"Wl^{x)Wl^{x"  -  Ax) 

'<5  { K? + <"  -  *">  ('  -  0]  +  (‘  -  0] } 

X  -  2  |l4(z')  [aPx^  +  (a;  -  x")  fl  -  0J  +  Vy{z')  Apyj  +  (p  -  y")  f  1  “  0  } 
-lA6k^  p  dz'C^(z')  J  d^f'  j  d^x"' Wilis' )Wl^{x"' -  Ax) 

^  ('  -  01’  ^  (‘  -  01T” 

X  -  2  {v^.(z')  [Ap,i + (X'  -  x”)  (i  -  0]  +  n(/)  [ap,£  +  iv  -  r)  (i  -  0]} 

+2.92k^  dz'C^(z')  J  dV  J  d^x"  Wl^(x')Wl^ix"  -  Ax) 
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-  2  |v;(2')  +  (x'  -  x")  (^1  -  +  Vy{2')  Apy^  +  {y'  -  y")  (l  -  } 


(3.100) 


Notice  that 


Vr(z') 


Ap.^  +  (x'  -  x")  (l  -  +  Vy{z')  Lpy^  +  {y'  -  y")  (l-pj  = 

V(z')-\Ap-  +  ix-x')(l-  -) 

L  2,  V  z,J 


(3.101) 


So,  equation  3.100  may  be  rewritten  as 


aC,(Ax,r), 
dr 


|t=0 


2'!  -i 


2\  -6 


2.43*2  y'*  Az'Cliz')  J  dV  J  dV'W^0(x)Vy0(x"- Ax) 
x{[Ap.i  +  (x-x»)(l-0]  +[ap,^  +  (p-!,")(i-0]  I 

xK(y)-  [Ap|  +  (x--x")  (l-0] 

+2.43*2  jT''  d2'C2(/)  J  d'x'  J  dV"  W^0(f')M^0(x"' -  Ax) 

X  {  [ap0  +  (X'  -  X"')  (l  -  0]  ’  0AP0  +  -  y'")  (l  -  0; 

xV(x-).[Ap0  +  (f-?”)(l-0] 

-4.86*2  y"'  dz'C2(x')  J  ^2^  J  d^x"Wl^{x')Wl^{x"  -  Ax) 

xV(z')  ■  Apj  +  -  x")  (l  -  (3.102) 


To  match  the  form  of  equation  3.82,  rewrite  equation  3.102  a.s 


gC,(Ax,  r) 
dr 


r=0  =  /  Az' Cl{z')V{z')-wv{z')  (3.103) 

Jo 
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where  wv{z')  is  the  wind  path  weighting  function  given  by 

Wv(z')  ^  l^2A3k^  '  dz'C^iz')  j  d^x  j  dV' -  Af) 

X  I  [Ap.i  +  (.  -  X")  (l  -  I)]  +  +  (y  -  f)  (l  -  £ 

+2A3k- dz' C^iz')  J  d^f'  J  d^x"'Wl^{x')Wl^ix"'-Ax) 

X  {  [ap£  +  (X'  -  X”)  (l  -  0]%  [ap£  +  (l  -  ^ 

-d.Sek^  dz' Cliz')  J  d^f'  J  d^x"Wl^{x')W^x*' -  Ax) 


y  \  1  2'1  3 


’'i  I  K'l + <*'  ■  *">  ('  ■  I)]  +  Kl + -  *'">  (‘  ■  1)1 }  1 " 


(3.104) 


where  f  is  a  unit  vector  along  Apj-  +  (?—  ^')  given  by 

[Apxff- Ax(l- [Ap,f^- Ay(l-ff)jy 
(l  -  f:)]'+  [APvff  -  Ay  (l  - 


(3.105) 


3.6.6  Reduction  of  V  Path  Weighting  Function,  wv{z')  The  steps  taken  to  simplify 
equation  3.104  are  identical  to  the  steps  taken  in  equations  3.31  through  3.44  of  the 
path  weighting  function  derivation.  The  result  is 

Wv{z')  =  i^2.43k^Li^  J  dw  tri 

+K(0  +  “(‘-0]  } 

’  { h"'  (y  -  ('  - 1)  -'''']+  h»  (0 + “  -  0]  1 1 


+2.43k^L2^  J  du  tri  ^ 


u  -  Ay" 
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"In  3  K  (0 + “  -  9. 

-  {h-  (0  -  ('  -  0  +  "■’’]  +  K  (0  +  “  -  0. 

-  I  h-  (0  -  ('  -  0  -  '^’'1  +  K  (0  +  “  ('  -  0. 


T  2  ^  3 


n  2  ^  3 


-A.Sek^L^^L^^ 

■  £i^+A„ 

dii  ^11  -|- 

£l^+Ay 

+  ^1 

du  5103(11)  -t-  ^ 

L  . 

JLX^  +  Ay 

i/2  +  i/J 


-  Ay^  gwaiu) 


u  +  Ay]gwz{u) 


(3.106) 


where 


gw^iu) 


■Ax  + 


-Ax  — 


,  ^”4. + (- 

-{h00 

z'  ( 

Apx—  +  l^-Aa:  + 


I^) 

Li  +  i/2 


-.2 

z'  (  z'\ 

+ 

K;;+“('-0J 

12^  3 


/  \  1 2'j  i 


+ 


2 

Ll  +  i/2 
2 

Li  —  L2 


)('-a]’-Kf“(-a]’) 

)('-;)r*hs-('-£)r) 


(3.107) 


Recall  that  both  Ay  and  Apy  were  set  equal  to  zero  for  the  path  weighting 
function.  This  resulted  in  no  separation  of  the  optical  beams  in  the  y  direction.  However, 
such  a  separation  is  crucial  for  measuring  wind  velocity.  The  underlying  phenomenon 
being  exploited  is  the  transport  of  eddies  from  one  ray  path  to  the  other.  If  the  ray  paths 
intersect,  it  is  impossible  to  infer  anything  about  the  rate  of  this  eddy  transport.  Let  the 
separation  in  the  x  direction  equal  zero  at  the  crossing  altitude,  therefore  some  offset  of 
the  ray  paths  in  the  y  direction  is  required.  As  shown  in  Figure  3.6,  this  offset  is  achieved 
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where 


X  <  2 


1_ 

L2 


-{ 


1 2  'I  3 


-i+Ay 


.(ii+ii)4.Ay 


(■ 


du  {  M  + 


.ii^+Ay  fl 

+  (-T- 


-  u  +  Ay 


>  r 

(3.108) 


gwziu) 


=  {Kf:+( 


Ax  + 

Ax  - 


L2  ~  Lj 

2 

Li  +  L2' 


)(-0]00'-a] 


/  \  T  2  ^  3 


-{K*{— i^X'-er-K'-aiT 

*  (Kr  •  (-“■  •  (‘  -  31'  •  K‘  -  91 T 


To  make  the  integrals  over  u  dimensionless,  make  the  changes  of  variables  Li  =  RL2  and 


u'  =  The  result  is 

Lq 


Wv(z')  =  |2.43fc^i2  ^L2^  J  dw'  tri  ^w'  —  ^ 

’  oftr-iX'  biT 


X  <  2 


-{I 


7  z'  { Ax  \  ( 


l2'\  3 


RL 


1 


-{iffi-*9ffi-)r*K-9r 
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4-2.43PZ2  ^  J  dw'  tri 


.....  7 

X  <  2  ^  1—  +  (  —  )  — 
La  V2j/  L2  j 


+ 


K‘-0] 


■.2'i  3 


i 

1 

2 

+ 

r 

[7  ,2'  fAx  V 

2 

'  ,  /  z'\ 

fl 

0-  -1-  1  +  -  —  -  1 

+ 

u'il - 

1 

.L2  V  L2  J . 

L  V  zj 

J 

-•1  r  /•^ir-+^ 

-4Mk'^R-n^  •  ^ 


d«' 


+ 


1  +  i?  Ay^ 


+i? 


/J 


■+f? 


J  /  /  /.  /-R  +  l  ,  AyX  ,  ,, 

dM  gw,{u  )  +  +  i;  j  ^ 


(3.110) 


where 


gwsiu') 


-{ 


\j_  ,  2'  C^x  l-R 

2 


2_  __  + 1 

1L2  2 


*f:(S-¥)l-|-'('-r)l 
•;ffi*^)i'*i'(-ai'l 


,2'l  3 


2^  3 


(3.111) 


where  7  =  Ap^rf — ^2;.  Equation  3.110  is  the  most  general  form  of  the  wind  path  weighting 
function.  However,  it  can  be  simplified  by  assuming  1  —  w  1.  The  simplified  form  is 


wv(z')  =  i^2A3k^R  ^L2^  J  du'  tri  ^w'  — 


X  <  2 


RL 


\  2 

:) 


ML. 


)2  ®  /  \  2 


+2.43it^Z,2  ^  j  d«'  tri  ^u'  - 


3-41 


X  <  2 


(i)  +“i  -[(i-')  -[(i+o 

)  gwsiu') 


/  du'  (u'  +  - 

y_F«±iU^  V 


-4MPR-^L2^* 

-+f? 


+  R  Aj; 


r^+ff  /R-\-l  Aw\  11 

Ji^+%^  V  2  Is/  JJ 


(3.112) 


where 


gw3(u‘ 


-1 


7  I  —  R 


T  2 


I2  ^  2 


•f-l 


7  ^  +  11 


+  u' 


-I 


r  7  ft  + 1 


t2 


Z/2  ^  2 


+ 


Again  recall  that  the  dependence  of  wv(z'}  upon  z'  is  contained  in  7. 


3.7  Form  of  the  V  Path  Weighting  Function 

The  key  result  in  extending  the  new  measurement  technique  to  the  remote  sensing 
of  transverse  atmospheric  winds  is  the  form  of  the  wind  path  weighting  function,  wy. 
Normalized  plots  of  the  x-  and  j/-directed  components  of  wv{z')  are  shown  in  Figures  3.7 
and  3.8.  In  both  cases,  the  value  used  for  the  aperture  separation  in  the  y  direction  was 
At/  =  y/L2-  Such  a  separation  could  be  easily  achieved  and  maintained  in  most  optical 
systems.  Figure  3.9  demonstrates  that  a  larger  separation  results  in  a  slower  decay  rate 
for  the  path  weighting  function. 

These  figures  illustrate  the  important  point  that  this  measurement  technique  mea¬ 
sures  the  wind  velocity  in  the  y  direction  only.  As  seen  in  Figure  3.7,  the  i-directed  path 
weighting  function  has  odd  symmetry  about  the  intersection  point  (7  =  0).  This  behavior 
occurs  because  eddies  above  the  intersection  point  shown  on  Figure  3.6  and  moving  in 
the  X  direction  intersect  the  beam  from  p2  first,  while  eddies  below  the  intersection  point 
intercept  the  beam  from  pi  first. 

The  j/-directed  wind  path  weighting  functions  shown  in  Figures  3.8  and  3.9  are  sym- 
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Figure  3.7.  i-Directed  Component  of  the  Wind  Path  Weighting  Function  wy 


y/i^f 

Figure  3.8.  j/-Directed  Component  of  the  Wind  Path  Weighting  Function  Wy 
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Y/Li 

Figure  3.9.  y-Directcd  Component  of  the  Wind  Path  Weighting  Function  Wy  for  'I'wo 
Values  of  A?/ 

metrical  about  ^  =  0  and  exhibit  rapid  decay  away  from  the  intersection  point.  Uoth 
pro|)erties  are  desirable  in  order  to  measure  the  y  com|)onent  of  the  wind.  Notice  that, 
just  as  for  the  6';’  [)ath  weighting  function,  the  horizontal  axis  of  the  plot  is  normalized  by 
Ij\.  Once  again,  the  ratio  R  ~  ^  plays  an  important  role  in  determining  the  shape  of  th(' 
path  weighting  function.  As  R  approaches  one,  the  wind  path  weighting  function  becomes 
more  sharply  peaked. 

.?.8  Rctiolnlion  of  V  Ralh  Weighting  Funclion 

As  was  done  for  the  path  weighting  function,  the  RMS  width  of  the  V  path 
weighting  function  was  computed  and  then  used  to  determine  the  path  weighting  function’s 
vertical  resolution.  Figure  .3.10  presents  the  re.solution  of  the  wind  weighting  function  as  a. 
function  of  the  source  .separation.  The  values  used  in  generating  the  plot  were  Ay  =  \/7^, 
R  =  0.9,  Ij\  =  10  cm,  and  a  calculated  RMS  width  of  the  path  weighting  function  of 
-y'  =  0.04.')  meters.  Reading  from  the  figure,  10  centimeter  vertical  resolution  is  possible 
using  guide  stars  .separated  by  4a  degrees.  It  should  be  noted,  however,  that  t  hese  results 
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6  (degrees) 

Figure  3.10.  Vertical  Resolution  of  V  Path  Weighting  Function 

are  based  upon  the  assumption  that  the  vertical  profile  of  is  perfectly  known.  If 
can  only  be  measured  to  within  some  resolution  Az',  and  its  value  is  changing  over  Az', 
then  the  resolution  of  the  wind  path  weighting  function  will  be  reduced. 

3. 9  Conclusion 

This  chapter  has  presented  the  derivation  of  the  and  V  path  weighting  functions 
for  a  new  crossed  beam  measurement  technique  which  relies  upon  correlations  of  band 
pass  filtered  outputs  of  wave  front  slope  sensors.  The  shapes  of  these  functions  have 
been  presented  and  their  theoretical  vertical  resolutions  calculated.  The  sampling  function 
nature  of  the  weighting  functions  indicates  that  this  technique  allows  vertical  profiling  of 
and  V,  with  vertical  resolution  dependent  upon  the  angular  separation  of  the  reference 
sources.  The  next  chapter  will  examine  the  ability  of  the  proposed  measurement  technique 
to  accurately  estimate  and  V . 
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IV.  Signal-to- Noise  Ratio  Analysis 


4-1  Introduction 

The  signal-to-noise  ratio  of  the  measurement  correlations  gives  insight  to  the  useful¬ 
ness  of  the  measurement  techniques  presented  in  Chapter  3.  The  quantities  of  interest, 
and  V,  are  derived  from  correlations  of  wave  front  phase  measurements,  (C,).  However,  C, 
is  a  random  process,  and  the  correlation  value  obtained  by  a  single  measurement  represents 
only  one  of  the  many  possible  values  of  the  correlation  function.  The  ratio  of  the  mean 
of  C,  to  the  standard  deviation  of  C,  is  a  good  indication  of  the  accuracy  of  the  estimate 
of  the  correlation  from  a  single  measurement  (24).  This  chapter  presents  the  derivation  of 
the  signal-to-noise  ratio  expressions  for  the  estimates  of  (C,)  and  (^((0)).  Additionally, 
signal-to-noise  ratio  results  are  shown  for  both  cases,  and  the  effects  of  measurement  noise 
on  the  SNR  are  examined. 

4.2  Derivation  of  SNR  Expression  for  (C,) 

4.2.1  Definition  of  {€,)  SNR  The  SNR  of  {€,)  is  defined  as 

SNR  = - - r  (4.1) 

where  KC,)!  is  the  magnitude  of  the  mean  value  of  C,,  and  (C^)  is  the  mean  square  value 
ofCV 

The  mean  has  been  derived  in  Chapter  3  and  is  given  in  equation  3.2,  where  Wc{z') 
is  found  in  equation  3.69.  The  remaining  task  is  to  find  an  expression  for  (Cf ). 

4.2.2  Derivation  of  {CD  The  mean  square  value  of  C,  is  defined  by 

{C])  =  ([(Sl  -  5,.  )(S2  -  S2'  )f )  (4.2) 

where  s„  is  the  output  of  slope  sensor  n.  The  expanded  version  of  equation  4.2  is 
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(4.3) 


(C;)  = 

((6162)^)  +  ((Sl'-S2')^)  +  ((^r'S2)^)  +  ((•51'S2')^) 

+  4(Si'S2'51'S2)  ~  2(.Si525i52')  ~  2(Si52Si'S2) 

—  2(5^  52' 5^  #52' )  —  2(5i'525i'52') 

As  can  be  seen,  the  mean  square  value  of  the  measurement  correlations  is  composed  of 
many  terms.  Results  for  each  term  will  be  presented  separately  and  later  combined.  Many 
of  the  assumptions  and  results  found  in  Chapter  3  will  be  referred  to  here.  Making  use  of 
the  slope  sensor  output  model  given  in  equation  3.10  and  beginning  with  the  first  term  in 
equation  4.3 

((5152)^)  = 

IJII  d2£dVdV'dV'W£,(x)VF^^(x'  -  Ax)Wl^{x'')Wl^{x'"  -  Ax) 

+  11  d^xd-VW/,(x)W7,(f")(<^,(x)^i(x"))(a202) 

+ 1  ydWx"'1^2,(x'-Ax)M^^.(x'"-Ax)((^2(f')<^2(n)(«i«i) 

+  (Q:iQ:i)(Q!2<12)  (4.4) 

where  a„  is  the  measurement  noise  of  sensor  n. 

Most  of  the  terms  in  equation  4.3  have  a  form  similar  to  equation  4.4.  In  the  interest 
of  brevity,  the  first  term  will  be  derived  in  some  detail,  but  only  the  final  expressions  for 
the  other  terms  will  be  presented. 

The  first  component  of  the  first  term  contains  the  fourth  order  moment  of  the  wave 
front  phase.  Assuming  Gaussian  statistics  for  the  wave  front  phase,  the  fourth  order 
moment  may  be  written  as  a  combination  of  second  order  moments  as  follows  (6) 

+{4>AS)M^")){MnMn)  (4.5) 
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Use  of  equation  4.5  in  equation  4.4  results  in 


((^1^2)  )  — 


I J  J  J -  Ax) 

X  [{(t>i{S)M^)){M^')M^"))  +  {4>i{^)M^')){<t>2{x')M^")) 


+  {<t>l{x)<l>2{^")){U^)<l>l{S"))] 

+  JJ  d=*f'dV"VU^^(f'  -  Ax)WlX^"  -  Ax){US^)<t>2(r')){a,a,) 
+(®i®i)(o;2<5;2)  (4-6) 


The  derivation  of  the  first  term  will  begin  with  the  evaluation  of  the  three  components 
arising  from  the  fourth  order  moment  of  the  phase.  The  first  of  these  components  is 


Comp  1  = 

I J  d^xd'^x'Wl^ix)Wl^{x'  -  Ax){M^)U^)) 
xj  y’dV'dV"VU/^(r')VU2,(f'"- Af)(</.i(x")<^2(^"'))  (4.7) 


This  component  consists  of  two  terms,  both  of  which  are  identical  to  equation  3.17.  Thus, 
the  expression  for  this  component  can  be  found  by  using  the  result  of  the  evaluation  of 
equation  3.17,  which  is  shown  in  equation  3.24.  Using  this  result,  the  first  component  is 


Comp  1  = 


d2:'C^(z')  Wi{z') 


2 


(4.8) 


where  Wi(z')  is  given  in  equation  3.47.  The  second  component  arising  from  the  fourth 
order  moment  of  the  phase  is 


Comp  2  = 

j  J  d^M^x''w[^{x)wiin{4>i(S)M^')) 

xJ  I  dWf"'VU/^(f'- Ai)VU£,(x"'- Af)(02(x')<^2(x"'))  (4.9) 
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Tho  terms  in  this  component  consist  of  correlations  of  the  wave  front  phase  from  a  single 
source  illuminating  a  single  aperture.  In  this  case  both  Ap^.  and  Apy  equal  zero  (there  is 
only  one  source),  and  the  aperture  separation  also  equals  zero  (there  is  only  one  aperture), 
so  the  appropriate  expression  is 


Comp  2  =  /  Az'  C‘^{z')  wa{z') 

Uo 


(4.10) 


where 


Wa(z')  =  —2.92pLi  ^  /  du'tri(M') 


Finally,  the  third  component  is 


Comp  3  = 

J  J  d'fdV'W^,(£)l^£,(x'"-A£)(<^i(x)<^2(f'")) 

xjj  dVd'f'W^^(x'-Ax)VF’,(x")(<^2(f'Vi(f"))  (4.12) 

This  component  is  identical  to  the  first  component,  the  only  difference  being  the  dummy 
variables  of  integration.  Thus,  its  result  is  also  given  by  equation  4.8. 

Moving  on  to  the  terms  containing  noise  in  equation  4.6,  the  first  noise  term  is  given 
by 

Noise  Term  ^  =  j  j  d"idV'VF£,(x)VF’^(x"){</>i(^)</>i(^'))(«2a2)  (4.13) 

Comparison  of  equation  4.13  with  equation  4.9  shows  that  this  term  can  be  expressed  as 

J  J  d^xd'^i?'Wl^{x)Wl^{x"){4>i{x)<f>i{x")){a2a2)  =  (^1 dz'  Cl{z')  Wa{z')  (4.14) 

where  (t|  =  (0302)  as  given  by  equation  3.14.  Similarly,  the  second  term  containing  a  noise 
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factor  is 


<Tf  f  dz'Cl{z')WA{z') 
Jo 


(4.15) 

Finally,  the  first  term  in  equation  4.3  may  now  be  written  by  combining  equations  4.8, 
4.10,  4.14,  and  4.15 


Cl{z')w,iz') 


+ 


dz'  C„(z')wa(z') 


dz'C^;(z')wA(z') 

r  dz'Cliz')wAz') 
Jo 


1  2  2 
+  (TICT2 


(4.16) 


where  Wa{z')  is  given  in  equation  4.11. 

Before  results  are  given  for  the  remaining  terms  in  equation  4.3,  there  are  two  ex¬ 
pressions  which  need  to  be  presented.  These  expressions  occur  often  in  the  remaining 
terms  but  have  not  been  encountered  before  in  this  thesis.  Therefore  they  are  evaluated 
now.  The  first  expression  is  the  ensemble  average  of  the  slope  measurement  noise  from 
two  co-located  slope  sensors  of  different  sizes,  and  is  defined  as 

(Q„a„.>  =  R^al,  (4.17) 

where  fi  =  a„  is  the  slope  measurement  noise  of  sensor  n,  sensor  n  is  a  square  with 
dimension  Li,  and  sensor  n'  has  dimension  L2-  The  variance  of  the  slope  measurement  noise 
from  sensor  n'  is  denoted  as  As  mentioned  before,  Wallner  has  shown  the  correlation  of 
noise  processes  from  separate  sensors  to  be  zero.  However,  in  the  case  present  in  equation 
4.17  the  two  subapertures  actually  overlap.  In  this  case,  the  correlation  yields  some  non¬ 
zero  value  which  is  dependent  upon  the  amount  of  overlap.  In  practice,  the  two  co-located 
apertures  will  each  be  composed  of  several  subapertures.  Because  they  are  co-located,  the 
two  apertures  will  share  some  number  of  subapertures.  The  more  subapertures  they  share, 
the  more  the  apertures’  noise  processes  will  be  correlated.  Equation  4.17  reflects  this  fact 
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through  the  factor  R-. 

The  second  term  not  encountered  before  is  the  correlation  of  the  wave  front  phase  of 
a  single  source  illuminating  two  co-located  apertures  of  different  sizes.  The  expression  for 
this  case  is 

J  j  d^(f)  d^(f")  dz'  Cl{z')ws{z')  (4.18) 

where 


wb{z') 


-2.92k^R~~L.  ^ 


J'L  ■*“'  (“' + “  fS 


/ 


-  “')  3bW) 


(4.19) 


where 


3b(w')  = 

R 


(4.20) 


All  results  necessary  to  evaluate  equation  4.3  are  now  available.  Expressions  for 
each  term  are  presented  and  then  combined  according  to  equation  4.3.  The  first  term  of 
equation  4.3  has  already  been  derived  and  is  shown  in  equation  4.16.  The  second  term  is 
given  by 


((Sl-S2')^)  — 


dz'  Cl{z')w2iz')  + 

r  dz'Cl(z')wA.iz') 
.Jo 

2 

r  dz'ci{z')wMz') 

Jo 

+  a^,  f  dz'  C^(z')wA'iz') 

.Jo 

1  -,2  _2 
“r(7|  /(T2/ 


(4.21) 
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where  11^2(2')  is  given  in  equation  3.50  and 
=  —2.92k^L2^  j  du'  tri{«') 


(4.22) 


The  third  term  is 


((-Sl"S2)^}  - 

2[|'’  +  d*'C'(z>^.(z')  !'■ 

+  <72  /  dz'Cl{z')WA-{z')  +  (Tl  f  dz' CI{z')wa{z') 
Jo  Jo 


dz"Cliz")WAizn 


(4.23) 


where  Waiz')  is  defined  in  equation  3.63.  Next,  the  fourth  term  is 


{{SlS^Y)  = 

2  [£'  iz'  +  j''  iz'CHz')w^,{z')  £ 

+0-2,  I  dz'  CI{z')wa{z')  d-  crl  f  dz'  C^{z')wA'iz') 

Jo  Jo 


dz"Cl{z")wA{z") 


(4.24) 


The  fifth  term  is 


{Si'S2'SiS2}  — 

dz' Cl{z')wi(z')  dz"  C^{z")w2{z")  +  dz' CI{z')wb{z') 
d^'C^(z')— j  R^icrl  +  al)  dz' Cl{z')wBiz') 


+ 

Uo 

+ftV^4 


(4.25) 
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where  Wb{z')  is  given  in  equation  4.19.  The  sixth  term  in  equation  4.3  is 

2^ )  “ 

2  r  dz' Cl{z')w,iz')  r  dz"Cl{z")^p- 
Jo  Jo  I 

+  r  dz'Cl{z')wB{z')  r  dz"c:,(z")w^iz") 

Jo  Jo 

d-R'^al  r  dz'Cl{z')w^iz')-\-(7l  r  dz' CI{z')wb{z') 

Jo  Jo 

+R''(Ti<tI  (4.26) 

Term  seven  is 

(S1S2S1/S2)  = 

2  r  dz' C^„(z')wi(z')  r  dz"C^(z")^fl 
Jo  Jo  2 

+  dz'Cl{z')wB{z')£  dz"Cl{z")w^{z") 

+al  r  dz'Cl{z‘)wB{z')-\-R^cTl  r  dz' Cl{z')w^{z') 

Jo  Jo 

^R?cyl(Tl  (4.27) 

The  eighth  term  is 

(5iS2'5i<S2')  = 

2  r  dz'C^Jz')^^  r  dz"Cliz")w2(z") 

Jo  A  Jo 

+  r  dz'Cl{z')wB{z')  r  dz"Cl{z")w^.{z") 

Jo  Jo 

+  <^2'  /  ^z:' C^{z')Wb(z')  +  R^<7i>  /  dz'  C^{z')Wai{z') 

Jo  Jo 

+  R^alal  (4.28) 

Finally,  the  ninth  term  is 

(Si<S2'Si'S2)  = 

2  r  dz'Cliz')"^  r  dz"Cl{z")w,(z") 

Jo  I  Jo 
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+  r  dz'C^^{z')WB{z')  r  dz"Cl{z")w^.{z") 

Jo  Jo 

+R^(tI  I  dz' C^{z')wa.{z')  +  crl  I  dz' C^(z')wb{z') 

Jo  Jo 


(4.29) 


Combining  the  nine  terms  in  equations  4.16  through  4.29  in  accordance  with  equation 
4.3  yields 

{CD  = 

2  dz'  Cl{z')w,{z')^  +  2  dz'  CDz')w,iz')\^'  +  dz'Cl{z')w^(z’)^' 

+  y"  dz'Cl{z')wAz')\  +  4  yy  dz'  CDz')ws{z')^\  8  d^'  CUz')"^^ 

+2  r  dz'Cliz')w^{z')  r  dz"CDzDwAiz") 

Jo  Jo 

+4  f  dz'  CDz')wiiz')  I  dz"  Cl(z")w2(z") 

Jo  Jo 

-8  /"  d^'cj(z')t«s(/)  r 

Jo  Jo  I 

-4  r  dz' CDz')WBiz')  r  dz"CDz")WA(z") 

Jo  Jo 

-4  r  dz' CDz')wb(z')  r  dz"CDz")wAiz") 

Jo  Jo 

+  y^  dz'CDz')wAiz')  4-  ^  dz'C^(/)«;^,(2r')J  [ct^  +  +  (af,  +  7|,)(1  -  2i?^)] 

+  [4fi^(cri,  +  cTj,)  -  2(<Ti  +(7l  +  CTj,  +  a|,)]  j  dz'  CI{z')wb(z') 

+  <Tf,<T|,(l  +  4R*  -  4RD  +  (<^^<7-2  +  -  2R^)  (4.30) 


At  this  point,  let  w'^  represent  w„  with  no  coefficients  before  the  integral.  For  exam¬ 
ple,  let 


iv'a{z)  = 


WAiz) 


(4.31) 


-2.92k^Li  ^ 

Making  such  substitutions  allows  the  coefficients  to  be  shown  in  equation  4.30  without 
having  to  show  all  the  accompanying  integrals.  Additionally,  let  Li  =  RL2,  and  the  result 
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IS 


(cn  = 


4.26k'^R-Hp  y^'  dz'  Cl{z')w[{z')  +  4.26k^LP  dz'  C^(z')w'^{z') 

+8.53k''R-iLp  yy  dz'C^(z')w'^iz')  +  8.53kHp  [jj"'  dz'Cl{z')w'^,{z') 

+34.11it^J^-'‘Z2^  yy  dz'  CI{z')w'b{z')^  +  n.05k^R-^L~^  dz' C^(z')w'3{z') 

4-l7.05k^R---Lf^  r  dz'C^(z')w'^{z')  f’  dz"Cl{z")w'^,{z") 

Jo  Jo 

4-8.53k^R-H;^  r  dz' C^(z')w',(z')  dz"  Cliz")w'^iz") 

Jo  Jo 

-limk^R-ii;^  r  dz' C^(z')w',{z')  r  dz"  C^{z")w'^{z") 

Jo  Jo 

-17.05k^R-H;^  r  dz'  Cl{z')w'^{z')  r  dz"  C^(z")w's(z") 

Jo  Jo 

-34.nk^R-H;^^  r  dz' Cl(z')w'giz')  r  dz"C^(z")w'^iz") 

Jo  Jo 

-34.11ib'/?-2£2^  I  '  dz' Cliz')w'B{z')  r  dz"Cl{z")w'^,{z") 

Jo  Jo 

+  \-2.92k^R~H;^^  r  dz'Cl{z')w'j,{z')-2meLt^  dz'Cl{z')w'^\z') 

Jo  Jo 

X  [trf  +  (tI  +  ((Tj,  +  Cj,)!  !  -  2R^)] 

-  [4R^{al,  +  al,)  -  2{a^  +  cTj  +  <Tj,  +  232k'^R~^Ll^  f  dz'  Cl{z')w'g{z') 

Jo 


(4.32) 


4.2.3  Reduction  of  Mean  Squared  Expression  The  resolution  of  an  optical  system 
operating  in  the  presence  of  atmospheric  turbulence  increases  with  aperture  size  until  a 
limiting  value  of  resolution  is  reached  beyond  which  aperture  size  increases  have  no  effect. 
This  limiting  value  is  called  the  atmospheric  coherence  diameter,  r-Q,  and  is  considered 
a  measure  of  the  quality  of  “seeing”  (6).  It  is  convenient  to  relate  the  SNR  to  Tq.  The 


expression  for  Tq  is 


rn  =  0.185 


/  C^„{z')dz' 


(4.33) 
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where  \  This  allows  fc'*  to  be  written  as 

k 


5.612 


[/  C2(^')d^f  " 

Use  of  equation  4.34  in  equation  4.32  results  in 


(4.34) 


^  _ A 

roJ  [/Q(. 


1  2 


.'L2 

+23.91  — 


+47.87  I  — 
+47.87 


roJ  [fC^ni 


+95.68 


VoJ  UC^^{z')dz 

L'“ 

■•oy  [/CJ( 
_ i;' 

)  (/cj( 

•^2 


/  r  \  3  T  r  *1  ^ 


u 

(Toot'*'*!’ 


-191.43  1 


+ 


-6.92 


To;  /Q(^')dr 

i  r  -2 


-6.92 


X2\5 


(t\  +  CTj^  +  ({T^/  +  CTjOCl  ~  2/?^) 
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(4.35) 


—  +  (Xy)  ~  2(cr(  +<72+  (7[i  +  <72<)j 


+  <7^(72^  +  <7(V2^(1  +  4R^  —  4R^)  +  (<7'i"<T2  +  <r'j"<72/)(l  -  2ii^) 


where  <  =  <7^ 

The  mean  of  C,  can  also  be  written  in  terms  of  Tq  using  equation  3.69 

(C.)  =  -3.46  (+)  ‘  J.  £  <i-''Cj(4')  [«-*»;(-'')  +  «-;(-'■)  -  2fi- vyz'i 

(4.36) 

Notice  that  the  factor  ^  J  common  to  equation  4.36  and  the  square  root  of 

equation  4.35.  Also,  recall  that  equation  4.36  is  in  the  numerator  of  the  SNR  expression, 
while  the  square  root  of  equation  4.35  is  in  the  denominator.  Therefore,  the  common  terms 
will  cancel  in  the  SNR  expression. 

All  that  remains  is  to  relate  the  slope  measurement  error,  (7„,  to  the  slope  sensor 
photon  count,  N .  This  will  allow  an  examination  of  the  light  intensity  necessary  for  useful 
results.  Welsh  and  Gardner  give  an  expression  for  the  slope  measurement  error  as  (23) 


0.867r77 

N^Tq 


L>  To 


(4.37) 


where  q  is  a  parameter  which  accounts  for  detector  imperfections  and  N  is  the  total 
subaperture  photon  count. 

Equation  4.37  allows  <7„  to  be  relaterl  to  <7„/.  Recall  that  C7„.  is  the  slope  measurement 
error  for  an  aperture  of  length  L2  and  <7„  is  the  slope  meas'.rement  error  for  an  aperture 
of  length  Li-  If  the  photon  count  for  an  aperture  of  dimension  L2  is  N ,  then  the  photon 
count,  N' ,  for  a  smaller,  co-located  aperture  of  dimension  L,  is 


r  2 

N'  ^  =  R^N 

L2 


(4.38) 


(T„  =  R  V„. 


(4.39) 
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Recall  that  a'^,  =  a'^,  Ll. 

It  is  possible  to  relate  a'^,  to  the  photon  noise  count  through  equation  4.37.  The 
resulting  equation  is 


(4.40) 


Apertures  .Sj/  and  S2'  a^re  identical,  therefore  it  is  reasonable  to  assume  that  a'^,  =  a'^,. 
Likewise,  it  is  reasonable  to  assume  that  o-f  =  a'^ .  Making  these  assumptions  and  using 
equations  4.38  and  4.40,  the  mean  square  value  of  C,  is 


23.917)!-^  dz' Cf,{z')w[iz')  +23.91  dz' C^iz')w'^{z') 

+47.87R-5  dz'Cl{z')w\iz')  +  11.95  dz'Cl(z')w'^,(z') 

+  l9lA3R-'^y^’  dz' Cl{z')w'Biz')^  +  95.68i?-''  J^'  dz' Cl{z')w's{z') 

+95.68i?-^  f  d/C^;{z')w'^iz')  r  dz"Cliz")wUz") 

Jo  Jo 

+47.87fl-5  r  dz' Cl{z')w[{z')  r  dz"  Cl(z'')w'^iz") 

Jo  Jo 

-95.68R-5  dz' Cliz')w[{z')  f’  dz"  C^(z")w'3iz") 

Jo  Jo 

-95.68R-'  r  dz' C^(z'}w',iz')  f'  dz" 

JO  Jo 

-191.43R-^  r  dz' Cl{z')w's{z')  r  dz"C^Jz")w'^iz") 

Jo  Jo 

-191.43R-^  r  dz' Cl{z')w'^{z')  r  dz"Cl{z"}w'^,iz") 

Jo  Jo 

(  iv  )"  f 

-13.84(0.867r77)2  - r  (T  -  2i?^  +  /  C^(z')dz' 

v(^)  / 

X  fi?-^  r  dz'Cl{z')w'^{z')+  r  dz'Cl{z')xv'^,(z') 

Jo  Jo 


-27.68(0.86x7/)'“*  ^  ,  [  Cl{z')dz'R~^(2R‘^  -  R-*  -  1)  /  ’  dz' C;(z')u4(z') 

l(^)V 
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+ 


(0.86177])^ 


/  C^Jz')dz' 

1(^)7  ^ 


(47?''  -  47?-  -  47?  +  27?-'  +  T?-  ’  +  1 ) 


(4.41) 


4-2.4  Simplification  of  SNR  Expression  The  SNR  expression  can  be  simplified  by 
making  the  assumption  that  1  —  ss  1.  Such  an  assumption  simplifies  both  the  mean  and 
the  mean  square  expressions.  The  expression  for  the  mean  square  value  of  C,  becomes 


(c;) 


23.91i?' 


f 

Jo 


dz'Cl(z')w';iz') 


+47.87(7?-=  +  1)+''  /  ‘  dz'Cl(z') 
Uo 

+ 191.43 /?-''fi2  r*  d2'C^(z')l 

Jo 


+  23.91 
2 


n  2 


+  95.687?-'' 

J 

+95.687?-5/12  ’  dz'C'2(2')] 

+47.877?-^  r  dz' Cl{z')w'l{z')  f'  dz"  Cl{z")w'i{z") 
Jo  Jo 

-95.687?-^  r  dz' C^iz')w'^iz')  dz"  C^{z")wU2'’) 
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where 


The  two  constants,  A  and  B,  found  in  equation  4.42  are  given  by 


and 


where  expressions  for  w'l{z'),  w^iyz')  and  w'^iz')  are  given  in  equations  4.43,  4.44,  and  4.45 
above. 

The  signal-to-noise  ratio  may  now  be  evaluated  using  equations  4.42  and  4.50  in 
equation  4.1. 


4-3  Signal-to- Noise  Ratio  Results 

The  signal-to-noise  ratio  gives  an  indication  of  the  ability  of  the  measurement  meth¬ 
ods  presented  in  this  thesis  to  provide  useful  estimates  of  C^.  In  order  to  evaluate  the 
SNR  expression  derived  in  the  previous  section,  it  was  necessary  to  choose  values  for  sev¬ 
eral  parameters.  After  choosing  the  necessary  parameters,  the  SNR  without  any  slope 
measurement  noise  present  (i.e.  N  oo)  was  calculated  for  various  ratios  of  Next, 
the  total  photon  count,  N,  was  allowed  to  vary  in  order  to  observe  the  effect  of  slope 
measurement  noise  on  the  SNR. 

The  model  used  for  the  C^iz')  profile  was  the  Hufnagel-Valley  model  shown  in  Chap¬ 
ter  2.  The  guide  stars  were  placed  at  z,  —  100  km  and  the  vertical  profile  was  only 
measured  to  an  altitude  of  z'  =  15  km.  This  allowed  use  of  the  approximation  1  —  ^  ss  1. 
The  guide  star  separation  was  Ap^  =  I  km  (corresponding  to  an  angular  separation  of 
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0.001  radians)  and  the  aperture  separation  was  allowed  to  vary  in  order  to  control  the 
intersection  altitude. 

4.3.1  SNR  Without  Measurement  Error  Figures  4.1  through  4.3  show  the  calcu¬ 
lated  signal-to-noise  ratios  when  no  slope  measurement  noise  is  present.  The  SNR  im- 


Figure  4.1.  SNR  Without  Measurement  Noise,  R  =  0.25 


proves  as  R  approaches  one.  This  behavior  indicates  that  the  variance  of  the  measurement 
correlations  is  approaching  zero  faster  than  the  mean  approaches  zero.  Another  way  of 
looking  at  it  is  that  the  weighting  function  is  becoming  more  sharply  peaked  as  ii  -+  1,  thus 
correlations  due  to  turbulent  eddies  not  at  the  intersection  altitude  (noise)  are  receiving 
less  weight. 

4.3.2  Measurement  Noise  Effects  on  SNR  As  equation  4.37  shows,  the  slope  mea¬ 
surement  error  is  related  to  the  photon  count.  Figure  4.4  demonstrates  the  effect  of  changes 
in  total  photon  count,  A,  on  the  SNR.  In  obtaining  this  data,  a  value  of  t)  =  1.5  was  as¬ 
sumed.  This  is  a  typical  value  of  tj  for  wave  front  phase  slope  sensors  (23).  The  y  axis  of 
the  plot  is  normalized  by  the  SNR  value  when  no  slope  measurement  noise  is  present,  i.e. 
the  SNR  for  A  -+  00.  The  figure  shows  that  the  lower  ratio  of  ^  allows  near  noiseless 
performance  at  a  lower  photon  count  value  than  required  by  the  higher  ratio.  This  is  due 
to  less  correlation  of  the  two  co-located  apertures’  measurement  errors  for  smaller  values  of 
Note  that  while  Figure  4.4  presents  the  SNR  as  a  function  of  noise  for  one  particular 
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Figure  4.2.  SNR  Without  Measurement  Noise,  R  =  0.5 


Figure  4.3.  SNR  Without  Measurement  Noise,  R  =  0.9 
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Figure  4.4.  SNR  as  a  Function  of  Noise 


altitude,  the  shape  of  the  curves  would  be  the  same  for  any  altitude  because  the  slope 
measurement  noise  is  independent  of  altitude. 


4-4  Derivation  of  SNR  Expression  for  (C'(0)) 

4.4- i  Definition  of  (C'(0))  SNR  The  wind  velocity  Y  is  estimated  from  the  slope 
evaluated  at  zero  time  lag  of  the  time-lagged  measurement  correlations.  Therefore  the 
SNR  of  the  estimate  of  C'  is  found  using  the  mean  and  variance  of  C'  at  r  =  0.  The 
definition  of  the  SNR  of  the  slope  measurement  at  r  =  0  is 


SNR 


lim 

=  T  — 0 


^jC.(r)-C.(0)'  j 


C.(r)-C.(0)j^^  _  |•^C.(r)-C.(0)^j^j  " 


(4.51) 


or 


_ (C,(r))  -  (C.(0)) _ 

{(Q(^)>  +  (^2(0))  -  2(C2(r)C2(0))  -  [(C,(r))  -  (C.(0))]2}  = 


(4. .52) 
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Expressions  for  the  mean  and  mean  square  values  of  the  measurement  correlations 
evaluated  at  r  =  0  have  already  been  derived  in  the  first  part  of  this  chapter.  However, 
equation  4.52  also  contains  the  mean  and  mean  square  values  of  C,  as  a  function  of  r, 
as  well  as  the  time-lagged  correlation  of  the  measurements.  Expressions  for  these  terms 
must  be  developed.  However,  the  derivation  of  these  remaining  terms  is  almost  identical 
to  previous  derivations  and  thus  is  not  presented  here.  Instead,  only  the  final  results  are 
shown. 

4-4-^  Mean  Value  as  a  Function  of  t  The  expression  for  the  mean  value  of  as 
a  function  of  r  is 
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where  6  is  the  angle  between  V  and  the  vector  f  defined  by  f  =  +  (z  -  x')  ^1  - 

and 
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(4.54) 

where  7  =  Apj.—  —  Aa;. 

2  s 

4-4-3  Expression  for  Mean  Square  Value  Terms  The  denominator  of  equation  4.52 
contains  the  sum  of  and  (C^(0)}.  {Cf(0))  has  already  been  calculated,  and  adding 

(CKt)}  yields 
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and 
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WU7^(2')  = 

''“'  (“'  +  ^  77)  ^  /jO0’ 

+  C(^-“'  +  !7)^-'"'><+-) 


where  fwz{u')  is  defined  in  equation  4.54.  The  weighting  functions  w[{z'),  w^iz'),  w'^[z') 
and  w'g{z')  are  defined  as  they  were  in  section  4.2.2.  Also, 

»7  =  /_'  du'  tri(.')  1 1  (1  -  0]  I  -  I  [,-£]+[„•(,-  0]  I  I  (4,59) 


As  always,  7  =  ^Prj-  -  Aa;.  The  only  expression  remaining  to  be  evaluated  in  equation 
4.52  is  the  time-lagged  correlation. 
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4-4-4  Expression  for  Time-lagged  Term  Before  the  SNR  can  be  evaluated,  an  ex¬ 
pression  is  necessary  for  the  time- lagged  correlation  contained  in  equation  4.52.  The  ex¬ 
pression  is 
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where  w\{z'),  w'^^z'),  w'^lz'),  ww[{z'),  ww^iz'),  ww'^(z'),  w\(z')  and  w'^(z')  have  all  been 
dedned  previously.  The  other  weighting  functions  present  in  equation  4.60  are  defined  as 
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where 


Coefficients  common  to  equations  4.53,  4.55,  and  4.60  are  not  shown  because  the 
common  coefficients  will  cancel  when  these  equations  are  placed  in  the  signal-to-noise 
ratio  expression  given  by  equation  4.52.  It  is  now  possible  to  calculate  SNR  values  using 
equations  4.53,  4.55,  and  4.60.  Before  doing  so,  however,  these  equations  are  simplified  by 
assuming  1  -  «  1. 

4.4-5  Simplification  of  SNR  Terms  Before  presenting  SNR  results,  the  terms  which 
comprise  the  SNR  expression  are  simplified  by  assuming  1  —  ^  «  1.  As  mentioned  in  Chap¬ 
ter  3,  such  an  approximation  is  often  valid  and  reduces  the  complexity  of  the  calculations 
which  must  be  performed  to  obtain  estimates  of  and  V .  The  resulting  simplified  ex¬ 
pression  for  the  mean  value  a.s  a  function  of  r  is 


{C,{t))  = -3.46  j  Az' Cf^{z')  R  ^ww"(z')  +  ww'fiz')  -  2w'W3(z')j  (4.65) 


The  simplified  form  of  the  sum  of  the  mean  square  value  terms  is 
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The  simplified  form  of  the  time-lagged  correlation  is 
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where  w'l{z'),  w'.^{z'),  and  w'^{z')  are  defined  in  equations  4.43,  4.44,  and  4.45.  The  con¬ 
stants  A  and  B  are  defined  in  equations  4.47  and  4.48.  The  other  weighting  functions 
present  in  equations  4.65,  4.66,  and  4.67  are  defined  as 
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and 
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Finally,  the  constants  M,  N,  and  P  are  given  by 
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Reiterating,  0  is  the  angle  between  V  and  Ap^  +  x  -  f'. 
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4-5  (<^'(0))  Signal-to- Noise  Ratio  Results 


The  previous  section  has  presented  the  expressions  needed  to  evaluate  the  SNR  of 
the  wind  sensing  technique  proposed  in  this  thesis.  In  order  to  evaluate  the  SNR,  values 
had  to  be  selected  for  several  parameters.  The  Bufton  wind  profile  presented  in  Chapter 
2  was  used  for  V{z'),  and  the  Hufnagel- Valley  turbulence  model  was  used  for  The 

assumed  measurement  geometry  was  the  same  as  when  evaluating  the  C,  SNR.  The  guide 
stars  were  placed  at  z,  =  100  km  and  the  V  vertical  profile  was  only  measured  to  z'  =  15 
km.  This  allowed  use  of  the  approximation  1  —  f-  ~  1.  The  guide  star  separation  was 
Apx  =  1  km  and  the  aperture  separation  was  allowed  to  vary  in  order  to  control  the 
intersection  altitude.  The  limit  as  r  — >  0  of  equation  4.52  was  then  evaluated  numerically. 

4.5.1  {C'(0))  SNR  Without  Measurement  Error  Figures  4.5  and  4.6  show  the  cal¬ 

culated  signal-to-noise  ratios  when  no  slope  measurement  noise  is  present.  Like  the  SNR 


Figure  4.5.  (C'(0))  SNR  Without  Measurement  Noise,  R  =  0.25 


of  the  C,  estimates,  the  SNR  of  the  C'(0)  estimates  improves  as  R  approaches  one.  The 
results  in  Figures  4.5  and  4.6  were  computed  using  cos  6  =  1. 

4.5.2  Measurement  Noise  Effects  on  (C^O))  SNR  Figure  4.7  .shows  the  the  effect 
of  changes  in  total  photon  count,  N,  on  the  SNR.  The  value  selected  for  7?  was  t)  =  1.5. 
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Normalized  SNR 


0 . 0025 


Figure  4.6.  (<^'(0)}  SNR  Without  Measurement  Noise,  R  —  0.9 


Figure  4.7.  (C'(0))  SNR  as  a  Function  of  Measurement  Noise 
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The  plot  is  normalized  so  that  an  SNR  value  of  one  is  the  normalized  SNR  value  when  no 
slope  measurement  error  is  present.  These  results  indicate  that  as  the  ratio  of  ^  becomes 
smaller,  less  photons  are  needed  to  achieve  near  noiseless  performance.  This  trend  is  due 
to  less  correlation  between  the  noise  processes  of  the  co  located  apertures.  Like  Figure  4.4, 
Figure  4.7  presents  the  SNR  as  a  function  of  noise  for  one  particular  altitude,  but  the  shape 
of  the  curves  would  be  the  same  for  any  altitude  because  the  noise  value  is  independent  of 
altitude. 

4.6  Improving  SNR 

The  signal-to-noise  ratio  never  exceeded  one  in  any  of  the  results  presented  thus  far. 
However,  these  results  were  all  for  a  single  correlation  measurement  on  one  aperture.  As 
mentioned  before,  slope  sensors  normally  consist  of  multiple  subapertures.  Thus,  it  may 
be  possible  to  make  more  than  one  measurement  at  a  time.  Also,  several  independent 
measurements  could  be  taken  on  the  same  aperture  over  a  period  of  time.  Either  method 
of  increasing  the  number  of  independent  estimates  available  will  increase  the  SNR  by 
decreasing  the  variance  of  the  estimate. 

4.7  Conclusion 

This  chapter  presented  the  SNR  expressions  and  results  for  the  estimates  of  the 
measurement  correlations  used  to  derive  and  V.  These  expressions  are  sufficiently 
general  that  any  desired  model  for  the  vertical  profiles  of  and  V  may  be  included. 
These  SNR  expressions  provide  important  information  for  the  determination  of  the  number 
of  independent  measurements  necessary  to  obtain  accurate  estimates  of  the  quantities  of 
interest.  Additionally,  the  slope  measurement  noise  was  related  to  total  photon  count  and 
the  effects  of  photon  count  on  the  SNR  quantified. 
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V.  Implementation  Issues 


5. 1  Introduction 

The  results  shown  in  previous  chapters  have  been  obtained  neglecting  certain  effects 
present  in  a  practical  implementation  of  the  proposed  measurement  technique.  As  a  result, 
the  performance  of  a  remote  sensing  system  based  upon  the  method  examined  in  this  thesis 
will  likely  be  less  than  the  theoretical  performance  presented  in  Chapters  III  and  IV. 

This  chapter  discusses  the  effects  some  of  the  implementation  issues  may  have  on 
system  performance.  The  chapter  begins  with  a  presentation  of  weighting  functions  calcu¬ 
lated  for  a  representative  measurement  geometry.  The  weighting  functions  are  calculated 
using  the  results  of  Chapter  3.  There  is  then  a  discussion  of  some  effects  which  will  al¬ 
ter  the  shapes  of  the  weighting  functions  and  hence  alter  the  performance  of  the  system. 
Possible  methods  for  dealing  with  these  effects  are  also  discussed. 

5.2  Example  Measurement  Geometry 

A  practical  measurement  geometry  is  presented  in  Figure  5.1.  In  this  example,  the 
guide  stars  are  placed  at  an  altitude  of  10  -  20  km  and  are  separated  by  1  -  3  m.  In 
this  case,  the  approximation  1  —  ^  w  1  is  not  valid.  The  aperture  in  the  sensor  plane 
is  assumed  to  be  either  a  1.5  m  x  1.5  m  aperture  constructed  of  16  subapertures,  or  a 
3.5  m  X  3.5  m  aperture  constructed  of  32  subapertures.  The  subapertures  are  assumed  to 
be  square. 

Figure  5.2  displays  the  path  weighting  function,  Wc{z'),  resulting  from  three 
different  sensor  separations,  in  this  case  z,  =  10  km,  Apx  =  1  m,  and  Lj  =  18.8  cm,  where 
Li  is  found  from  Li  =  2  Figure  5.3  shows  the  j/  directed  wind  weighting  function 

Wy(z')  for  the  same  conditions.  The  figures  clearly  demonstrate  that  the  resolution  of 
the  weighting  functions  decreases  at  lower  altitudes.  Recall  from  Section  3.5  that  the 
resolution  is  inversely  proportional  to  the  angle  at  which  the  optical  beams  intersect.  This 
angle  decreases  as  the  intersection  altitude  decreases.  Increasing  the  altitude  of  the  guide 
stars  would  help  mitigate  this  effect. 


5-1 


Altitude  z 


Sensor  Plane 

Figure  5.1.  A  Representative  Measurement  Geometry 


Figures  5.4  and  5.5  present  the  weighting  functions  when  z,  =  20  km,  Ap^  =  3 
m,  and  Li  =  2  (^  |~)  =  22  cm.  Once  again  the  resolution  of  the  weighting  functions 
decreases  with  decreasing  intersection  altitude. 

The  weighting  functions  in  Figures  5.2  through  5.5  were  calculated  assuming  sta¬ 
tionary  guide  star  images  in  the  image  plane  and  deterministic  guide  star  locations  within 
the  atmosphere.  In  a  real  implementation  of  this  measurement  technique,  neither  of  these 
assumptions  is  true.  Thus,  the  effects  of  image  dance  and  beam  wander  must  be  discussed. 
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Normalized  ‘*y(Z)  Normalized  n^(2i) 


Normalized  Normalized  ^(7-) 
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Figure  5.4.  Second  Example  of  wciz')  for  Three  Values  of  Aa- 


Figure  5.5.  Second  Example  of  wv{z')  ■  y  for  Three  Values  of  Ax 
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5.3  Tilt  Correction 

The  turbulence  induces  random  changes  in  the  overall  wave  front  phase  slope,  or  tilt, 
which  causes  the  guide  star  images  to  move,  or  dance,  about  the  image  plane.  Removing 
the  low  order  distortions  of  the  wave  front  phase  yields  a  pseudo-stationary  image  in 
the  image  plane.  In  most  measurement  schemes  incorporating  artificial  guide  stars,  this 
compen'^ation  is  usually  accomplished  with  a  tip-tilt  correction  system.  The  system  is 
composed  of  a  wave  front  tilt  sensor,  a  tip-tilt  correction  mirror,  and  a  control  algorithm. 

The  band  pass  filtering  approach  adopted  in  this  thesis  does  not  require  such  a 
system,  however.  The  low  order  phase  distortions  are  common  to  both  co-located  sensors 
and  are  removed  when  the  larger  sensor’s  output  is  subtracted  from  the  smaller  sensor’s 
output.  Thus  the  low  frequency  components  of  the  phase  distortions  do  not  contrib'ite  to 
the  measurement  correlations. 

5.4  Beam  Wander 

The  artificial  guide  stars  are  generated  by  ground  based  lasers.  As  a  laser  beam 
propagates  through  the  atmosphere,  turbulent  eddies  which  are  large  with  respect  to  the 
beam’s  diameter,  D,  deflect  the  beam.  This  causes  random  movement  of  each  guide  star’s 
location.  Thus,  the  guide  star  separation,  Ap,  becomes  a  random  variable  and  therefore 
the  crossing  altitude  of  the  beams  also  becomes  a  random  variable. 

If  the  uplinks  are  separated  by  a  distance  much  greater  than  the  beam  diameter, 
then  the  movements  of  the  guide  stars  are  uncorrelated  because  the  beams  are  affected  bv 
unrelated  turbulent  eddies.  The  movements  may  be  correlated  at  some  time  lag  due  to 
motion  of  the  turbulent  eddies.  The  amount  of  correlation  would  depend  upon  the  time 
lag  and  the  vertical  wind  profile. 

Fante  (2)  shows  that  the  variance  of  the  guide  star  position  is  dependent  upon  the 
ratio  of  where  po  is  called  the  spherical  wave  coherence  length.  Beland  and  Krause- 
Postorff  (1)  show  that  =  ^,  where  Cq  is  defined  in  equation  4.33. 

To  examine  the  effect  of  beam  wander  on  the  example  measurement  geometry  of 
Figure  5.1,  let  Cq  =  10  cm,  which  represents  moderately  good  seeing  conditions,  and  let 
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the  laser  beam  have  a  diameter  of  5  cm  and  wavelength  of  A  =  0.5/xm.  Assuming  the 
guide  stars  are  generated  at  =  10  km  and  the  crossing  altitude  is  5  km,  application 
of  Fante’s  equations  gives  the  variance  of  each  guidestar’s  movement  as  0.5  cm,  and  the 
RMS  value  of  the  intersection  altitude  deviation  is  22  m.  A  larger  ratio  of  ^  would  yield  a 
greater  crossing  altitude  deviation.  Obviously,  this  movement  of  the  intersection  altitude 
will  reduce  the  vertical  resolution  by  flattening  and  broadening  the  peaks  of  the  weighting 
functions. 

A  method  for  reducing  beam  wander  is  to  use  an  adaptive  optics  system  to  shape 
the  wave  front  so  the  beam  focuses  to  the  same  point  in  space  at  all  times.  Such  an 
approach  requires  knowledge  of  the  quantit.j  being  measured,  therefore  the  variance  of  the 
intersection  will  only  be  reduced  to  some  limiting  value  beyond  which  further  reduction  is 
not  possible.  The  amount  of  reduction  in  the  variance  of  the  intersection  altitude  would 
depend  upon  the  capabilities  of  the  adaptive  optics  system  employed. 

5.  -5  Conclusion 

This  chapter  presented  an  example  measurement  geometry  and  showed  the  resulting 
weighting  functions  as  calculated  from  the  expressions  in  chapter  three.  This  chapter  also 
discussed  the  effects  of  guide  star  motion  upon  the  performance  of  a  practical  system 
implementing  the  proposed  measurement  technique.  Methods  which  may  mitigate  these 
effects  so  that  the  performance  approximates  that  shown  in  the  example  were  suggested. 
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VI.  Conclusions  and  Recommendations 


6. 1  Conclusions 

This  research  has  demonstrated  that  the  spatial  correlations  of  band  pass  filtered 
wave  front  slope  measurements  may  be  used  to  calculate  the  value  of  the  atmospheric 
structure  constant.  The  spatial  correlations  are  related  to  through  an  integral  expres¬ 
sion  containing  and  a  path  weighting  function.  An  expression  for  the  path  weighting 
function  was  developed,  and  it  was  found  that  the  shape  of  the  path  weighting  function 
is  determined  by  the  measurement  geometry  and  by  the  width  of  the  filter  pass  band.  As 
the  pass  band  becomes  more  narrow,  the  path  weighting  function  becomes  more  sharply 
peaked. 

This  research  has  also  shown  that  the  temporal  correlations  of  band  pass  filtered 
wave  front  slope  measurements  may  be  used  to  obtain  the  transverse  atmospheric  wind 
velocity  profile.  Evaluation  of  the  slope  of  the  temporal  correlation  function  at  zero  time 
lag  yields  an  integral  expression  containing  C^,  V,  and  a  wind  path  weighting  function.  An 
expression  for  the  wind  path  weighting  function  was  developed,  and  again  it  was  found  that 
the  shape  of  the  path  weighting  function  is  determined  by  the  measurement  geometry  and 
by  the  width  of  the  filter  pass  band,  with  the  wind  path  weighting  function  becoming  more 
sharply  peaked  as  the  filter  pass  band  becomes  more  narrow.  This  method  measures  the 
wind  in  one  direction  only,  therefore  the  measurement  geometry  must  be  properly  aligned 
to  obtain  the  wind  speed  in  the  direction  of  interest. 

The  measurement  technique  presented  in  this  thesis  has  not  been  experimentally 
implemented.  However,  signal-to-noise  ratio  calculations  for  commonly  accepted  and 
V  profile  models  indicate  that  multiple  measurements  will  be  necessary  to  obtain  accurate 
estimates  of  the  quantities  being  sensed.  The  number  of  measurements  required  will  likely 
not  be  so  great  as  to  render  this  method  impractical  for  measurement  of  C^{z).  However, 
low  signal-to-noise  ratios  for  the  estimate  of  the  slope  of  the  temporal  correlation  function 
at  zero  time  lag  indicates  that  this  technique  may  not  be  practical  for  the  remote  sensing 
of  transverse  atmospheric  wind  profiles. 
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6.2  Recommendations 


The  following  additional  research  is  suggested. 

1.  Using  accepted  models  for  atmospheric  turbulence  and  wind  velocities,  this  measurement 
technique  should  be  computer  simulated.  The  spatial  and  temporal  correlations  of  filtered 
wave  front  slope  measurements  should  be  calculated  by  the  simulation,  and  the  results 
used  to  determine  and  V .  The  values  of  C'^  and  V  found  in  this  manner  should  then 
be  compared  to  known  results  to  determine  the  accuracy  of  this  method. 

2.  The  shapes  of  the  path  weighting  functions  as  modified  due  to  laser  propagation  effects 
should  be  determined. 

3.  This  optical  remote  sensing  technique  should  be  experimentally  implemented  and  the 
results  compared  to  those  obtained  from  existing  methods. 

6.3  Summary 

The  spatial  correlations  of  band  pass  filtered  wave  front  slope  measurements  may 
be  used  to  calculate  the  value  of  the  atmospheric  structure  constant.  Additionally,  the 
temporal  correlations  of  band  pass  filtered  wave  front  slope  measurements  may  be  used 
to  obtain  the  transverse  atmospheric  wind  velocity  profile.  Experimental  implementation 
and  verification  of  this  technique  is  necessary  before  its  use  in  an  adaptive  optics  system. 
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